INTERMOLECULAR  POTENTIAL  ENERGY  FUNCTIONS  FOR 
PAIRS  OF  SIMPLE  POLYATOMIC  MOLECULES  AND  THEIR 
APPLICATION  TO  PROPERTIES  OF  GASES  AND  LIQUIDS 
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PHEFACE 


Tli«re  hBV«  b«eB  several  attempts  to  celcuUte  poteotial  eoergy 
fuoctlons  for  relatively  small  polystaaie  moleovles  in  order  to  deter- 
mine  tbe  proper  fora  of  these  poteatiels.  The  end  result  of  these  oal- 
oulations  has  usaally  heea  the  oorrelatioo  of  properties  of  these 
polyatoaio  molecules  usiog  the  eslouleted  form  of  the  poteotial  energy 
function  as  the  basis.  These  cnloulations,  therefore,  have  been  more 
goalitative  chan  quantitative  io  their  applioatloo. 

In  Che  present  vork,  these  calculated  potential  energy  functions 
are  tested  in  a more  quantitative  manner.  In  Part  I,  Internetion  peram- 
eters  for  constituent  attma  of  polyatouio  molecules  are  determined  from 
inert  gas  and  dlatomio  gas  data,  and  are  then  used  in  the  oaloulstlon 
of  potential  energy  functions  for  the  more  complex  polyatccic  molecules 
considered.  Theie  celculated  potentiel  energy  functions  are  then  tested 
by  calculating  second  virial  coafficioats  and  oomparing  with  «perimental 


Host  of  the  results  of  cslculated  potentiol  eoergy  fuootions  i 
previous  Horks  hsve  boon  applied  to  tbe  gas  state.  In  Part  n of  I 
prosnot  work,  an  attempt  is  made  to  determine  if  the  aalculated  pot 
tlel  energy  functions  cf  Part  I con  be  nppiled  to  the  pr^erties  ni 
pure  liquids.  Beosuae  of  tbe  emplexity  of  using  olaberete  liquid 
modela.  the  present  work  makes  use  of  a simplified  model  of  the  lie 
state,  the  smoothed  potential  noil  model-  Experimental  mi 


tha  tharmsl  pressure  poefficieat  for  turee  ccapounds  are  preseatod,  aad 
the  molar  voliae,  themal  pressure  coefficient  uod  configurational 
energy-volma  relationsbips  ace  discussed  ia  light  of  the  predictloas 
of  this  model.  Also  lacludod  in  Part  II  is  an  investigution  of  the 
oonflpuratlcael  energj-voluaa  ralstlonsbip  in  liquids  using  the  radial 
ui  function  as  the  baoi 


Tba  author  ia  indebted  to  the  Natiooal  Science  foundetioo  for  pro- 
viding Cooperative  Graduate  Fellouahips  to  the  author  for  the  period 

been  mode.  The  author  also  thanhs  tbo  Dniveraity  of  Floride  Cenputlns 
Ceater  for  the  doautioa  of  coaputer  time  for  seae  of  the  calculations. 
The  author  appreolates  the  invaluable  help  of  Dr.  I.  U.  Heed.  1X1.  who 


directed  the  reaeerch  end  was  chairmen  of  the  supervisory  ccnmitteo. 
ond  the  oonilderetion  of  Dr.  I.  L.  Balloy,  m,  Dr.  A.  A.  Broyles, 

Dr-  I.  E.  Haxfield,  Dr.  J.  T.  Uoore.  snd  Professor  B.  D.  'dellcar.  fr.. 


in  the  supervisory  osaittee.  Thanks  are  also  di 
'e  for  her  potienoe  end  support  and  for  typing  tt 


Dyal  McKinley 


, CAICIWTIONS  FOS  POLTi>IO«C  KCmjJLES  . . 


V6 

CtCH3)4 


SI <063)4 
. SUMMRy  MO 


VI.  BASUEBENI  OF  (eP/8T)y  ill  ton  IMStWIDHES 


frOD  Csloul&led  PotontlBls  . . . . 
R /£  j:i*  r*/^  for  tb«  Spheriool 
Beduoftd  Boylo  Teaperature  for  tbe 


Shell  PotontlBl 
Spherlosl  Shell 


Potential  Enerjj  PiuiotiM  of  CjlCBjJj  

Pcteatlal  Energy  Function  of  SF^  

Potential  Energy  Punctiooe  of  SiF-  

Potential  Energy  Pnootion  of  Si«H3)4  - 1 

Tberael  Preeeure  Coefficient  and  Cnapreaeibility  of  CFa  . . . 

Ttemal  Presaure  Coeffloient  of  Ca^  

U^id  Molnr^olinte  Correinlona  and  Extrapoiated  Liipiid 


Configurational  Energiaa  and  tfolar  Voilaoea  of  Llgulds  . 

Paranetera  Uaed  in  Figures  43  and  44  

Belationsbip  Batuean  r*  nod  r"/^  

Spherical  SSell  Potential,  s 1.70  

Spherical  Shell  Potential,  r^  r 1-60  

Spherical  Shell  Potential,  r,  c 1.90  

Spherical  Shell  Potential,  r^  c S.OO  

Sphorloai  Shall  Potential,  rj  a 2.10  

Spherical  Shell  Potential,  rj  e 2.20  

Heduced  Boyle  Teaperaiuro  for  tie  Spnerlcal  Shall  Poieot 
Lannard-donea  6-12  Potential  . 

Lannard-Jonee  6-14  Potential  . 

Lannard-Jones  6-16  Potaotlal  , 

Lennard-Jonea  6-18  Potential  . 

Lannard-Jones  6-20  Potential  . 

Lennard-Jones  7-14  Poteotlal  . 

Lennard-Jonea  7-21  Potential  . 


Lsaaard-Jones  7-2fl  PotantlBl  

Thanaal  Preaaure  Coefflcleat  of  CCI^  

Tliormal  Praaaare  Coeffiolant  of  HeaaureaeDta  f end  2, . 
Taemal  Preaaure  Coefficient  of  CF>,  Vaaaureaenta  2 and  4, . 
Tberoal  Preaaure  Coaffieient  of  CFa,  Meaautementa  5 end  6.. 
Tberoal  Preaaure  Coefficient  of  CF,,  Neeaurenenta  7 and  8., 
Tbersal  Preaaure  Cceffioient  of  CF.,  tfeeauruientc  9 and  10. 
Tbereal  Preaaure  Coefficient  of  CF4,  Uesaurenenta  II  and  12 
Tlienaal  Preaaure  Coefficient  of  (7.,  Heaeurttaeota  13  and  14 
Tberfflal  Preaaure  Coefficieet  of  C8,,  Keaaurenenta  1 and  2. . 
Theraal  Preaaure  Coefficient  of  CH4,  Vaaeureaenta  3 end  4.. 
Tliemal  Preaaure  Coefficient  of  CH4,  Heaauranenta  S nod  6.. 
Tfieraal  Preaaure  Ceeffioient  of  CH4,  KeesuraieAtn  7 and  8,, 
Tbenaal  Preaaure  Coefficient  of  CH4,  ifeasuraneota  9 and  10, 


IlfTSODUCTIOK  AND  BB/IR 

is  aecssssrjr  to  have  soae  understaodiog  of  iatereotioas  os  the  nozeoular 
scale.  The  singlest  systm  io  Nhlcb  iaceraotioas  betweea  atoes  sad  doI- 
eoales  oafl  be  described  Is  cae  wblcb  oouists  of  an  isolated  pair  of 
atone  (or  BOleoalesl.  From  guantun  oieobenlcal  derivations,  it  has  bssa 
deternined  that  e pair  of  atoms  oil]  attract  eeoh  other  at  reUtivelF 
large  separatiCLOs  sad  nill  repel  each  other  at  very  saall  aeparatlons. 

At  scBO  intarTaediate  point  there  ivtu  be  no  net  force  of  iatersction. 
These  Interaotioos  are  described  by  the  nutual  pair  potential  eaergy  as 

live  of  the  slope  of  their  potential  eaergy  fuaetlon  at  that  seperetioo. 
The  exact  fom  of  tho  potential  eaergy  function  oannct,  in  general,  be 
derived  fran  quantum  r.eehanion  benausa  of  the  ccmtplexity  of  the  probien. 
For  this  reason  aeverol  sttenpts  hsve  been  mode  to  detennine  s nnthe- 
mntioaliy  simple  approxlDStioB  to  tbe  potential  energy  ft 
mould  at  least  approach  physical  reality. 

Poteotlal  Energy  Functions  for  Single  Centers  of  Interaction 
The  most  vridaly  used  of  these  approxlBstions  has  been  the 


Figure  1,  Lonnard-Jonau 


Potaotial 


free  roratioa  of  cbe  oolecules.  Tte  iotereotioos  between  ooosti'.Heat 
atosie  on  sdjaoeot  ooleouleo  »ero  asoiKoii  to  follow  tie  Unoirdnlooes 
6-12  potential  onergy  fonotion. 

The  ptocedBre  of  ooloolation  wat  aa  follows:  an  ewpresiioo  was 
written  for  the  enorgy  of  ioteraotloa  botwesn  every  pair  of  aicaa  tone 
fraj  each  otIbcuIo)  in  two  Interactlag  woleeules  as  a foaotSon  of  ori- 
entation and  distanoe  betMaea  nolecular  oeaters.  This  awpreasioa  was 
then  avoraged  over  all  orientations  of  the  two  ooleoulea,  giving  equal 
weight  to  all  orleotatioos.  The  total  mutual  nolaoular  pair  potential 
energy  was  oquol  to  the  sum  of  the  interaction  energies  for  the  eton 

The  integrations  in  the  averaging  pronass  ware  done  in  olosad  forwt 
whlnh  made  the  resnlta  easier  to  use  than  the  Infinite  series  solotion 
hy  Atoji  and  Upscowb  (4>  which  will  he  discussed  Ister. 

Althongh  the  results  were  epplioibie  to  any  two  moleoulas  of  the 
type  AjSji  and  0^0^,  Htiann  and  Laabort  made  ealoulationa  tor  ooly  ona 
oasa.  Tha  assunptioaa  used  wore  that 


batweeo  moleooles  of  the  aaae  hind; 


la  the  molecules  is  the  sage;  l.g. 


3.  the  wlniDton  in 


4.  r = , where  p is  the  distanoe  frwa  the 

cole  to  the  oeator  of  the  perlpharel  itows. 

Wheo  the  potentiel  energy  function  (a  function 
T-nl ecclar  oenters)  was  colculohed  for  this  cosa,  it 


Lwinoni-JonM  tjpa  of  poteoUol  onerjy  function  (Equation  2)  nlth 
exponents  m s 7 an«  o o 28,  JH„  a 7-28  potentSel  onnrBy  function-  Ilio 
nr  tois  potential  enarg?  funotlon  for 
(rial  coefficient  data.  Ttls  potential 

found  to  give  side  paranaters  for  ttese 
I froD  the  Lonnard-Jooea 


several  icoloeulea  fror,  neoi 
■as  found  to  correlate  gas 
and  C(CR3>^  vary  well,  and 


8-12  potential. 

given  bare.  The  follovdng  equations  apply  ti 
eoleculas  and  assionlng  only  that  t 
rotation  and  tbat  tbe  total  Interaction  ener( 
aotloos  between  atoq  pairs.  Equal  weight  la 
of  the  two  BOleculos.  No  other  aaauaptlous  » 

n(B)  = u^i-  U^B  + UBc  + ugj, 

where  D(R>  e Eotal  interaction  potential  oner 


seen  two  molecules  are 
ir  general  equatlona  will  be 
interactions  between  two 

tiveo  to  all  orientations 


B t Disteoce  batwoeo  uoleculer  centers 

e Potential  energy  of  Interactioc  between  central  scorns 
(atona  A and  C are  at  center  of  moleculeJ 
“ad*  % - Potential  energy  of  Interaoclon  between  central  and  all 
poripheral  atnss 

Uyi  = Potential  energy  of  interaction  betiioon  ell  peripheral 


(£) 


of  caoleeule  to  perlphoral 


p r DiitoBce  Iron  oenter 


T , £ n PsTSBietero  is  Leoaard'^cneB  6-12  potoatlel  for  atoas 
Subsoripta  A,  B,  C,  D rofar  to  t&o  atoas  iorolvaO. 

Notlea  that  theao  e^atlooa  oao  be  uaad  direetl;  far  iotaractlaoa 


OIB)  e lice  + + Ufp 


la  raloe  of  pcp  b«  “t®  the  C - F bead  length,  other  loleenles  of 
ita  typo  Inolude  SF;,,  SIF4,  CCMslo.  end  SKCBjl^,  If  CHj  groepa 
ra  ooneidared  to  ioteraot  as  single  oeotera. 

The  results  of  Haasnn  and  Laobert  can  alas  be  applied  to  ■oleoulas 
I Bhioh  there  ia  no  eton  at  the  center,  auoh  as  in  the  diataalo  gases 


Hj.  Ng.  Oj.  F2.  and  Clj.  For  these  noleoulas.  Equation  8 glues  the 
total  lutersotloo.  The  parar.etera  tor  these  molecules  are  a e o a 0, 
b e d c 2,  and  the  value  of  ^ in  each  ease  la  one-holf  the  bond  length 

There  la  another  type  of  noleoule  for  nhich  the  equations  of 
Hoiinnn  and  Uaiberl  can  be  uaad.  Those  are  svoleoulss  such  as  CjFe  and 
CjKBglj  fa  ehioh,  uhen  freely  rototing,  the  coastltnont  otoas  ace 


eoneentrlo  spheres.  Equetioa 


io  Intaraeciong  be 


is  datersiaed  tij 
iribuSB  BOtblag  l 


pair  poteotial  energy  funotlon  of  C2F'e-  * 


id  of  Equation  a only.  Eguntions  5,  6,  and  7 eon- 
tia  total  potantinl  energy  for  this  type  of  tioIocuU, 
tebieh  bos  no  atnn  st  Its  cantor.  If  It  is  desirable  to  caloulote  inter- 
aatloes  between  ■olecolee  In  wbieb  tlore  are  core  tbsn  two  spberiosl 
shells  of  atoms  uhea  the  Bolecule  Is  rotated  freely,  Egiiatioo  9 can  be 


tbe  interaoticas  aooi 
I esreful,  howeyer,  t3 


)f  tbe  BOleoQle  to 


to  atOBB  Involyed. 


lated  molecules  ec 


ItoLl  and  Ueiecrh  fi'i.  Atoji  and  LipsconA  also  assumed  that  tbe 
•a  of  istoraotiOB  between  two  polyatomic  moleeulea  are  In  the  con- 

;here  was  un-correlated  free  roietioa  of  all  uolocules,  and  that 
itel  mutual  potential  energy  of  iateraotion  betweeo  a pair  of  iao- 
>e  detonninad  as  the  sum  of  Interactions  of  all 
:iRS  avoragad  with  egual  weigbi  over  all  orienta- 
10  molaonlea.  In  order  to  robe  oslculetions.  tbe  authors 
le  oonstltnent  atoms  interact  according  to  the  lennard- 
Jones  6-12  potential  energy  funetioo.  Thus,  their  results  should  be 
equiveloct  to  those  of  Hanano  and  Lambert  since  tbe  eosumpclons  ere  the 
same-  This  wss  shown  to  he  true,  even  though  the  iategratlooe  wore 
performed  by  a differtnt  method  frcei  that  of  Ei 
Infinite  serine  solutlona. 


CalealscioQs  of  tbo  loiniBUm  lo  the  poteotisl 
CCI4,  and  SF^  eere  node  by  aoalgning  iotareotion 
e‘:llaent  atoao  and  uaing  the  derived  eguatlona. 


jy  fanciion  of  CF4, 


la  vlacoaity  ol 


energy  fuootloa  pi 


flnorine  atc=c  ia  CF4  aad  £ 
r i detemirted  frooi  date  01 
paraaeteri  for  oarboa  and  for  oulfar  atore,  alas. 

be  ohloriae  atons  vt 
! coaffioient  of  Clg  gaa. 

the  noleoalea  Bare  eoapared  with  6-12  potential 
re  detenlned  from  erparSisetital  data  aad  were 

al  energy  funotiona.  It  will  be  ahowa  later 
that  a cwgiBrlaon  of  £ aad  r'  valuta  wllhciut  taVlcg  the  abape  of  tie 
potential  energy  fuaetion  into  aoocunt  la  of  doubtful  va-ue. 

fla^pspu  Beleaoa  darlred  by  uae  of  Gegoabauer  polynoniale 

the  general  eapreaslon  for  tbe  interaotioa  of  an  leolated  pair  of 
spherloally  eyinetricti  polyatsilo  aicleoulee.  Re  nonaldered  tbe  oen- 
tare  of  Interection  to  be  in  Uie  oooatltoent  attraa,  and  alnce  tbe 
largeat  interaetiona  ere  between  paripberel  atiot,  be  conaidared  tbe 
entiro  ioteraotion  between  naleoolea  to  be  due  to  interootlona  mong 
peripheral  atnna.  Re  aaeuced  that  interaotioni  between  pairs  of  atoaa 
could  bo  ropreaonted  by  a Leonard-Jonao  6-12  potential  energy  funotloa 
and  that  tbe  total  Interaotion  between  uoleoulea  ie  the  ainn  of  the 
ioteraotions  between  all  pairs  of  atoma.  His  derivation  la  valid  for 
two  moleeulea  la  a given  position  and  doaa  not  atiaapt  to  give  tbe 
average  interaotion  over  nil. orientatlcat  of  tbe  two  noleculos  as  in  tbe 


iHo  disousB&d  ptevionsl;.  In  bis  gnnsrtL  expression  there  nre 

oentrel  terms  iteras  ebieb  do  not  depend  upon  Tooleeulsr  ibepe,  struc- 
ture, end  orientation)  and  ncn-oentral  terns  xhieb  depeod  upon  the 
shape,  structure,  and  crleotetion  of  the  moleenles,  Slece  those  non- 
central  terse  deiemine  the  crlentatlon  effects  for  polyetnsic  mole- 
cules, these  terms  mere  of  the  greotesc  Interest.  The  anther  coeoluded 
tbtt  for  syraetrScal  pclyntomie  moleeules  the  eon-centril  terms  mnhe 
only  a smell  contribution  ae  long  as  the  seperstioo  betueen  nclecules 
Is  large  compered  to  the  eiae  of  the  molecule.  The  central  terms  mere 
foond.  In  general,  to  agree  mlth  trcaimente  mhicb  assumed  free  roteticn 
of  the  molecoles  to  get  eversga  potnntlals.  Frsoi  the  date  on  the 
second  virUl  eoofflclont  of  and  the  central  toms  of  his  derivt- 
Cloo,  he  detarmlhed  that  the  Lenoard-Jonas  6-12  potential  Interaction 


f* 


For  CClg  be  also  fouad  that  £qi  ^ 

This  indleatei  that  lAea  ealeulatlog  interactions  bntxeen  perlpb- 
eral  at^s  in  molecules,  the  r*  values  for  the  ooble  gas  of  slmllor 
siao  can  probably  be  used  for  perlphorsl  atcas,  but  scmo  care  must  ba 
exerolsed  In  the  chclco  of  the  values  used. 


It  methods  foi 

deteminlng  potential  energy  feuiotlons  of  gases  la  the  scattarlog  o] 
Ians  la  gates.  Ibeao  expsrliuhts  usually  yield  tbs  potential  energi 
Interaotlon  over  very  norrcil  rsogos  of  eeporntions,  hoMaver.  For  tl 
reason  these  results,  at  the  prese.nt  time,  are  of  limited  value  in 


Seteiralcias  suffUiaat  portions  of  tin  potentlol  ooerjy  funetlons  to 
oslculBte  pb^sloal  propertloo.  It  Ic  hoped  that  this  will  not  bo  tree 
in  the  reor  future.  Tho  paper  by  Clcney  aad  Vauderslioe  covers  a 
narrow  range  of  ooleoular  separations,  but  sowo  of  their  oaUelotlons 
ore  simrlar  to  the  work  ooae  here  and  should  be  oentioned. 

Fro*  ioo  sootterlog  data,  potential  energy  funetlons  were  deter- 
ailned  for  the  sjstens  a*-  CB4,  H*-  CgSd,  [§  - CB4,  aj  - c^R^,  n+- 
®2^6'  ^ These  potential  energy  fonotjono 

ware  valid  nt  separations  whore  the  forces  were  repulsive.  The  authors 
used  the  potential  energy  functions  for  CS,  and  CP4  with  tie  hydrogen 
ions  to  predict  potontlel  energy  functions  for  hydrogen  ions  nod  C.K. 
end  CgPj. 

In  the  first  scheme,  ths  Interestlont  between  H*-  CjSj  and  a*-  C.F, 
Kara  assionad  to  be  due  entirely  to  Intoreetlons  between  hydrogen  Ions 
and  the  perlphnrsl  hydrogen  end  fluorine  atoms  U CjRj  a.d  CgF^. 
raspaoilvely.  The  Intarantlon  centers  in  the  peripheral  hydrogea  and 

Sana  assuaptidn  was  neda  for  Intereotlona  between  H*-  CR4  and  B*-  CP4. 
and  Interaction  poreteters  for  hydrogen  and  fluorine  stoos  with  B'*’were 
evaluated  froa  the  «■*■-  CH4  and  H*-  CF4  data.  Ihia  was  dona  by  deriving 
ecmations  for  tho  intaraoticn  of  eo  non  with  a freely  rotating  poly- 
atoalo  BOlaoula  as  in  the  papers  discussed  previously-  Those  poramaters 
for  hydrogen  and  fluorine  detarsiinad  Xrau  CH.  and  CF.  data  wore  tUcn 
used  to  oaloulote  tha  potential  energy  of  Interaotion  betwee.-.  CpH;, 
and  H*-  CjFj.  The  ealoulated  and  awporljiantal  results  sgreod  within  a 

oednre  was  followed  to  oaloulete  the  potential  energy  functions  for 


^ - CjHj  and  - Cgfj  tde  potantdal  anargy  funetiona  fi 


Tf.B  otlsr  whew  nwd  »as  to  ataieio  that  In  CjHj  and  CgFfc  t 
Earn  two  cecters  of  interaction.  Tor  esampla.  e< 

Paranatara  for  CEg  and  CP3  groups  oara  datamlned  as  before  from  date 
on  CHg  and  and  the  potential  onergy  functions  for  C2Pg  and  C2E^ 
acre  ealcuietcd.  The  agreeaer.t  again  was  within  a facicr  of  two 


cuAnai  II 


THE  SECOH)  nmtt  COEPnClENI  AKO  POTEMTIiU.  ENEfiGV  FDKTIONS 
The  equatloa  of  state  for  Inperfeot  pesos  oas  be  erittea  Id  the 


sa . 1 + 61U  + + esu  + . 


(W 


P z Absolute  pressure 
V a Volume  per  mole 


The  ooeffioiests  8<T>,  C(T).  and  DCT)  ere  oalled  the  seoond,  third,  end 
fourth  viriaL  coefficients,  reepeetiuely,  aad'are  dlreotly  related  to 
the  Interactions  between  the  molocnies  of  the  gas. 

This  equntloa  of  state  diverges  at  the  doneitj  of  the  llliultl  end  Is 
used  pritnarlly  for  gases  at  low  to  siediiun  deosltj,  Tho  coefficient  B(T) 
is  Introduoed  to  ocoonnt  for  interoctions  of  pairs  of  BOleoules;  C(I) 

four  moleoules,  etc.  Therefore,  at  low  density,  use  of  the  second 
virial  coeffiolant  is  usually  sufficient,  but  as  the  density  Inoreeses, 
other  virial  coefficients  must  also  be  used. 

There  Is  a large  volums  of  data  available  on  tbe  second  virial 
coeffioiant  of  gases,  bot  date  on  the  third  and  higher  virial  ooeffi- 
cients  are  not  as  plentiful  or  at  reliable  because  of 


experimental 


.]A. 


J3> 

«Se«b'''h,n)  US) 

'“T^  OD 


By  slBPl8  B6i;ipulatiui, 


.!=(!]■’  <13 

Even  for  sinpla  palyavanic  nolaoulas,  Cbe  ceabsrs  of  iateraotlcn 
Bill  differ  froo  one  another  both  in  looition  end  la  magaleude,  and  if 
the  potential  energy  funoiione  were  erpirlcally  fitted  by  tbe  generai- 
Ued  Leaaatd-Jonea  potautlel,  then  different  moleonlea  «ould  here 
different  relnea  of  d and  o.  It  is  aboun  in  Situation  14  that  the  aeoond 
ritlal  ooefflelent  is  a function  of  n and  n ns  oell  as  of  tempereture. 
Caleulacioae  of  the  saoond  virial  ooeffioient  oora  nada  by  hand  in  the 


sensitive  a funotlon  the  second  virial  ooeffioiant  ie  of  d and  n.  For 
onnvaniance  and  to  sehe  the  results  sore  useful,  tho  oaloulations  uere 
done  in  reduced  form.  The  reduced  second  virial  ooeffioiant  ia 


r-K 

torparature  T 1 


nhare  Ic  z Boltamann'e 
Caloulationa  of  E 
potectlal  energy  funot 


a fuootion  of  T*  have  haan  ni 


given  by  airrohfaldor,  Ciirtl», 


7-28  poten- 


Bird  {11 . 8nd  cdleulstioap  of  ^ es  b fiiBctlon  of  T*  for  the 
tSal  »sro  tiodo  bj  Hraojui  ooa  Ud>ort  (3).  Band  oaloulotlono  uoro  sodo 
Jn  tbis  proooM  MOTS  for  the  8-18,  6-24,  7-14.  oad  7-21  potantUl 
aaorgy  fanctlone  in  order  to  have  Doro  8 ii.  T*  oorveo  for  oBaporlson. 
Tbo  ioficlto  soriao  in  Eqaatlo,*:  14  caivorges  rapldlf  or  hlijh  voluoB  of 
l" , but  divotgo*  ot  Iw  T . Tho  esleuloted  B*  voluos  are  llmltail  for 
loa  raluea  of  T*  boOBOoe  of  tbo  ourier  of  terao  it  voo  noceoxor?  to  uxo 
to  got  reasonobU  aoouraoy  in  B*  for  loa  T*. 

Slooo  tile  geoerol  oqnationx  for  BQ)  for  an  m-n  potential  (Egua- 
tlraa  14-18)  era  rttber  oosplea.  Uio  final  equations  for  oalcuUting  b’ 
are  gluau  beloo  for  oaob  of  the  potential  onorg)'  fanotiona  BOntloned 


j = 4(2]-s/^  a*)'®''* 
T-ai  Pomillal 


'l»)-iZ[r(¥|]iTl 
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'’N-inrft 


Ihs  rasulu  of  tbeaa  hi 


id  caloalfftiaas  are  given  in  i 
7-40  pountuis  taken  freo  U 
to  BboGm  in  ?igore  2. 


The  aaoond  virlal  caeffioiant  ia  s 
form  of  the  potential  eoerg;  function  b 
cfteo  be  correlated  egualiy  noil  by  aev 
of  tie  paraxietare  4^  and  r®,  of  course 


lauee  anperiniental  data  can 
ral  w-a  potaotiels.  The  values 
depend  upon  the  potential  to 
is  of  interest  to  enaBine  the 
for  fariouc  ir,-n  potentials.  It 


ieh  the  data  ess  fitted.  Therefore, 
letionsbip  esong  the  £ and  ralu 
of  especial  interest  to  exsnine  the  vslues  since  they  ate  more 
ititlve  to  ohangoa  in  potential  enorgy  funetioo  than  r*  veluea. 

Tie  rolatlonehlp  among  c°  values  determined  from  arperlmental 

il  ocafficloni  deu  for  various  m-n  potentials  can  he  deter- 
' ctaalnstion  of  the  Boyle  point  for  m-o  potactlalt.  The  Boyle 


virial  ooaffieient  ii 
•Iso  aero.  From  the 


end  virial  aoeffioient  plota  given  li 
.-educed  Boyle  tenpareture  for  eaoh  ai 
ioce  for  the  potentials  shoiin  in 
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= 2.07  t..  r..u«^  .e«o.d  vi 


T«:per9tures  for 


Boduced 


hsviivg  &s  zairazt  value  c£  1.77,  Since  the  retie  of  ^’s  for  these  teo 
potentials  aheulfi  be  about  2.07,  »o  find  that  should  be 

Prou  this  examlaatioo  of  the  second  virlal  coefficleot,  as  oaloa* 
luted  from  several  difforeot  potentisl  energr  fonetlone,  one  oao  see 
that  altoough  the  data  aay  be  empirloaUp  fitted  to  several  potential 
aoorgy  fimctlons  equally  uell.  the  parameters  obtained  tfill  be  quite 
differeot  for  the  different  potentials.  Therefore,  one  should  exeroise 
oore  in  c.hooslne  the  potential  energy  funotion  to  ehlch  experinental 
second  viriel  coefficient  data  will  ba  fitted,  eapeoially  if  the  peram- 
etecs  are  to  he  used  to  celculate  ether  properties. 

Spheiiosl  Shell  Xedel 

One  of  the  limitations  upon  using  an  n-n  type  potential  energy 
function  to  fit  data  on  polyatonlo  gosaa  is  that,  In  geoerel,  such  a 
function  does  not  fit  e calculated  potantlel  energy  function  over  the 
antiro  ra.'.gQ  of  seperetions  of  interest,  and  tables  of  second  virial 
eoeffloiente  are  not  available  for  many  different  m-n  pcteotials.  Like' 
vise,  correlations  of  second  virial  ccefficieots  csing  n-o  type  poten- 
tial eaergy  funotions  are  not  easily  made  ever  vide  rsages  cf  tempera- 
ture with  the  seooed  virial  coefficient  ttblea  now  available. 

A recent  peper  by  da  hocoo  aod  Hoover  fd)  la  quite  useful  as  a 
step  in  removing  this  limitation.  They  hava  derived  equations  for  the 

culated  eutecaivs  tebles  by  the  use  of  an  electronic  computer  which 


UQifdrx  dIctributSoa  ci 


cadel  asea  Its  aas'^ptioiv  taat  tliere  ia  a 
LKinard-5ona*  6-12  lataraction  sites  aver  the 
spheres,  that  is,  bosause  of  free  retetihh  el 


•faee  of  two  ioteract'ar 
jpie  polyatceaie  nole- 
)n  the  surfeoe  of  a 


id  if  all  crfeotstions  are  ^Ivea  equal  neiqht,  the  total 
le  peripheral  atons  of  two  such  aoleoulas  is  given  b; 
the  interaction  of  the  too  spheriasl  strfeoes  which  have  a unfforv 
disttibutioh  of  Interaction  sites.  This  Intaraotion  describes  the 
spherical  shell  potential.  Notice  that  this  potential  describes  exaetli 


CPs,  since  all  of  the  atoos  DUSt  he  la  a aingla  spherical  ehell  for  this 


tl  to  opplf.  This  Is  eneotlp  the  loae  oese  es  covered  by  Equation  d, 
Ls.  Bith  thin  potential,  the  perar.ater  which  detennlnoa  the  nhape 


results,  Jis  r^/^  decreases,  the  well  in  the  potential  energy  funotioo 


ptraoeter  for  tho  shape  of  the  potential,  Xq.  The  paranstar  is 


defiaed  by  de  Boceo  and  Soever  es 


'■=2^ 

where  S is  the  separation  at  which  tht  Dininua  occurs  in  the  totsl 
potential  energy  for  the  ieolated  pair  of  noleculas.  Tha  authors  ware 
cainly  interested  in  better  oorrelation  of  expBrir.ental  second  virial 


ooefflcieot  data  end  used  do  electronic 


coloulcte  tho  host 


ir:s'er  sxMe  :<3turally  ia  thfi  derivation  of  the  oqtation 
id  viTial  ooeffioient  bj  da  Socco  and  Hoovar,  and  itoa  con- 
le  in  ibeir  caicrlationa  end  in  the  corrsietion  of  experi- 

eeioulate  the  seoond  virial  cooffloient  for  a given  mole- 
ita  the  apnerioal  aheU  model.  A more  conveaient  parameter 

and  ^ ia  the  diatance  from  the  tenter  of  the  molecule  to  the  peripheral 
atone.  Fortunatelv.  by  caiculatlng  several  apherltal  aheU  potentlala. 
one  tan  determine  the  relationahip  betneen  r*/^  aad  S /2^.  Several 
aphsritel  shell  pttentiale  for  different  values  of  r^  here  taltuleted 
in  the  present  horh,  end  froi  these  oaltuletions  the  relationship 
betmaon  R and  r i;aa  determined.  This  rolitionahip  is  given  in 
Table  A1  and  Figure  A1  of  Appendix  A. 


function  for  a molecule  in  Hhith  there  uori 
ptaltiooa  other  than  is  a single  spherical 
shell  potential  mould  net  fit  axectl:'.  In 
tally  idilch  spboricel  chell  potential  beat 

oaXculated  in  the  present 


culating  the  aphorloei  shell  poten- 

re  uore  interaotlon  eeotara  in 
.1  then  the  spherical 
:n  order  to  detartiine  eapiri- 
it  fits  the  oeloulated  potontial 
these  spherical  abell  pcttatial 
apherloal  shell  pctentlals  were 
init  tccputer  using  Eguation  S, 


and  are  given  in  tnbulsr  form  in  Appendix  A. 

The  reduced  Style  to^eraturo  es  a function  of  r^,  gi 
Table  AC  and  Figure  A2  of  ^pendix  A,  for  the  spherical  al 
‘.!aa  obtained  fran  the  second  virial  ocefficieots  given  ia 


It  S$  i2tK6rti.-.s  xe  sate  that  there  U very  little  aiKeraooe 
iettieee  the  T-28  potential  sod  the  iphevieal  shell  potential  tor 

the  tpheriool  shell  nnSel.  This  Is  shown  tj  the  oosporlson  of  seoend 
vlrisl  eoef'ioler.u  sheen  in  Figaro  4,  shore  SQ/SIrira'S,-!  j, 
against  T*  for  the  potentials  eentinaed  shove.  The  solid  lines  are  for 
the  sra  potentin’  and  the  lodioeted  points  ere  for  the  apherioal 


Fnplriosl  Egnetioos  for  the  Second  Vlrlel  Coefficient 

3T‘helpt.Emi8tiDH  fl.9).  The  Berthelot  egnetios  is  a cwo-parasi- 
eter  egnotloe  of  state  which  is  quite  useful  for  oslculeting  sesend 
vltiel  ooeffioients.  The  seooad  viriol  ooafficient  Ircj  this  egustion  is 


Pltner-I  Corr-TotUn  fid' 
perties  of  oonpolar  liquids  am 


ised  on  en  snpiricel  scentrio 
This  fsoter  ess  ohosen  to  take  ncccunt  of  the  deviations  oi 
il  energy  functions  of  reol  noleoulos  fro,i  the  h-l2  pctea- 
'Ideal  flnlo"  nhieh  would  fit  the  6-12  potentioi  nodal  has 


ooipleu  and  dcvi. 


potential, 


frcj  Tapor  pres- 


i.ura,  oritiial  tor^eratura,  an4  critical  prcacure  Jor  aay  nonpfliBr  ooa- 
pound.  Valnaa  of  tte  ecaotrio  faetor  oaleulatad  Of  Pltier  for  aOine 
oolaoules  and  iy  tuia  luthor  for  otfcers  ora  jlran  in  Taile  15,  elont 

Tea  aacnad  olrlal  coaffioienc  frai  taic  carralatlac  it 
= (0.1M5+  0.073a)  - (0.330  - 0.«glT“l 


- (0,1385 +O.SOU)t;2 

- o.ooTsur;:® 


Li^ccnb  (4],  wer«  derived  vhioh  eee  be  ueed  to  fialeulste 

poreotlal  eaergy  fodotiou  for  polyetodiia  noleaules.  Haneen  and  Lei- 
berb  and  de  Boccn  and  Bnover  used  rheee  reeultB  to  CQrrels&e  axpari*- 
naatai  eenood  virial  ooefflclecb  date. 


.aj  paremttara.  Ttay  did  net 


d ita  validity  ae 


potaatiala  bitb  exparlr.antsl  b-22  paten 

Xt  seems,  tberadare, 
fanotla.es  baa  not  been  fully 
beyond  qualitative  Investiqations. 

In  tbe  next  abapte.'S,  tbia  ealcule'.ioa  pr&Oudure  is  furtber 
7be  equetlons  derived  by  Bssaer  end  Lenbert  arc  used  to  calcelets 
pctentlal  energy  funotinns  Cfquetioae  d-Ii  of  tble  uorbl. 

Xn  tble  ohapter  tbe  potential  a.neruy  fanotiona  of  N^,  0^ 
and  CI2  are  oaXoalated  sod  eoi^ined  mitt 


7nese  persrteters  ere  nacded  for  oaicLlatione  for  loie  complex  molecules. 
Hote  Chet  these  molecules  fit  tbe  spberlcel  shell  rudel  exactly,  and 
that  Equation  8 can  bo  used  to  osloulsie  tbe  potential  energy  functioss. 


NotiBB  ttst  if  r*  is  know  for  th«  ooostituant  noenvs,  and  p U 
known  for  the  noleaulo,  then  froa  Equation  3.  one  c&n  oalculate  the 
ratio  3/f  without  koowlns  (3ecall  that  U Is  the  total  poteotlEl 
energy  of  lateraeiioOf  and  Is  the  eneroj  parmeter  for  the  oonetitu- 
snt  etonst)  The  cpjsntity  ^ is  sarely  one-half  the  bond  distSDoe 
between  stons  in  the  moleculB.  The  tobies  of  Bowen  (11)  contain  this 
information , and  all  needed  bond  lengthe  are  givsn  in  Tnble  4.  The 
site  peraoaters,  r*'s,  for  ths  atoms  in  these  moleoulea  should  not  be 
very  muoh  different  from  tboce  for  the  closest  ioert  gss.  Therefore, 

nitrogen,  ohygen,  and  flucaine  atons,  r”  will  be  taken  ns  that  for  neooj 
sad  for  ohlorlne  stns.  r*  will  be  taken  as  tbat  for  srgoo.  Porsrbtsrs 
for  tbs  inert  gases  for  the  b-12  potontial  are  given  in  Table  9,  along 


values  of  r and  tbe  ratio  con  be  calculated  fran  Equation  B 

as  a fonction  of  the  distance  betweeo  noieoular  oentere,  B.  When 
0/S  la  plotted  against  B,  one  con  looats  the  minimum  in  the  curve, 
which  will  heve  the  oonrdinalea  (B*,-E  /^*),  where  -E*  is  the  mlnimtfn 
value  nf  D.  Using  the  veluee  of  the  ooordiooies  at  the  minimion,  plots 
of  U/E*  :CM.~  E/B*  can  bo  prepared.  By  oc.oparing  these  reduced  potentiel 
energy  ourves  with  those  for  m-n  potentials  (given  In  Appendix  B),  ane 
can  find  an  m-o  potentiel  which  given  a ressooebla  fit  ti 


^6-12  ^ potentiel)  is  available  for  these  molecules 

from  expsrimeotal  gas  riocosities  and  sacood  vlrial  eoofflciencs,  end 
can  bo  used  to  determine  by  moltiplying  by  the  proper  footer  from 


Figaro  B- 


Figare  6.  Potenlial  Energy  Functions  of  Dlntemic 


flustin*  Odd  cttloiSne  alias  ware  ra-ealoulaiad  nsUg  ia«  spSarietl  Jliell 
goteolial.  Fr«  tie  apherloal  sieU  potmtlsl,  the  oaloulated  paraBewri 


A)j.  = M-7»iC 
(iS*/k)5j  I ail's; 


(O-T**  f«  lie 


the  la-D  potentiels  were 
fluorlee  atoa  and  aOS'K  tor  tbe  ohlsriae  alM.  Th) 
errer  ot  about  3 paroeut  In  the  energy  paraRetera, 
rcatiag  tie  calculated  apherloel  shell  poteatials  wj 
rtia  error  is  smaller  thee  tbe  uncertainty  in  the  I 
Hhieb  they  vrere  derived,  and  Here  not  large  enough  to  eerrant  repeating 
the  oaloulatieos.  Therefore,  later  calculatioca  in  which  these  paraiB- 


;b  ei-n  potentials. 


The  values  of  i*/k  tn  the  eonstitueat  atems  were  found  to  be  con- 
siderably larger  than  those  tor  the  corresponding  inert  gases  as  shoHh 
in  Table  7,  This  Is  in  contrediotlon  to  the  findings  of  Baleecu  <S). 

If  tbe  calculated  potential  energy  function  doss  not  differ  greatly 
fron  the  b-12  potential,  the  oalouloted  valhes  of  R should  be  close  to 
the  n”  values  for  tbs  b-12  potential  determined  from  experimental  data. 
This  can  be  seen  to  be  true  from  Table  7 for  the  diatonic  gases.  The 
Dost  way  to  test  the  values  of  It  calculated  for  the  diatofliic  molecules 
is  to  ealeuicte  second  vlrial  coeffloiants  using  these  valuas  of  S . 
Since  the  E values  used  in  the  above  calculetions  were  teben  frem 
experimental  data,  the  calculated  values  of  the  second  virial  ccoffi- 
ciont  will  be  correct  ot  the  Boyle  point.  If  ft  is  much  in  error,  how- 

corslderably  different  from  tbe  expetlmentsl  values  st  other  temperatures 


beoause  the  eeoond  Tlrial  eoeffloient  Is  prapartlonal  to  These 
ealculetlons  nere  nede  end  ere  described  seperetelp  beloe  for  each  of 


the  gesae.  Celcslaticns 


corrections  li 


Nitrogen  (Ngl 

coeffloieot  date  areilable  for  Ng.  EnperLeental  oeasumente  nere  feeds 
by  Uolborn  and  OUe  (19);  HioheU,  HoaterSt  end  do  Boer  C14);  Keyes  (15) 
Vhlte,  er  «1 . (16):  Reesoa  and  ven  Lemaerln  (17);  and  van  Icterbeek, 
et  al . (IB),  Tlieae  data  agree  eery  nell,  but  all  points  available  could 
not  be  need  because  of  their  nuaiber.  Enoug)i  points  were  used  to  shofe 
clearly  where  the  second  virlal  coefficient  carve  Is.  The  data  of 


Oxygen  (O2) 

en  »1.  (16>.  Tee  Boyle  tenperoturo  woo  delorainod  froa  nbe  data  of 
Byero  to  bo  403%.  Previously,  was  caleuUted  to  be  3.07,  aod  frcm 
Figure  A2  ttie  reduoed  Boyle  teaperatoro  is  2.64.  Free  tbose  values 
£*/]c  was  detemined  to  be  132.6%,  B*  wes  previously  oelculotod  to  be 

Figure  3.  The  egreaaent  caa  be  sees  to  be  BXoellent.  wbicb  indioatea 


Fluorlae  CF^l 

oxperlDontal  aaaaursients  of  the  second  virisi  ooefficlaat  of 
were  aade  by  HQlte,  et  al.  <16),  but  tbeir  values  are  probably  larger 


ia  tlisoluto  value  thau  tbe  actual  cascnil  vtrlal  ci 
tuperaturea  becauaa  tbey  tfld  net  use  tbe  tbird  virial  coefrieient  In 

ualuuUta  tba  eecond  virial  coefilclent  froa  tbe  6-ld  patential  tablee 
prscentM  earlier  and  frea  the  apherical  abell  poteutial  tablea  of 
da  Socco  and  Bouvsr  for  r^  = 2,74.  Tbeae  calculated  curves  are  con'* 
pared  In  Figure  4 with  the  experlnental  date  and  nlth  curves  colculated 
fr«  Pitret's  ccrrelatlon  and  the  Berthelot  equation,  ubicb  were  dis- 
cussed earlier.  It  can  be  seen  that  the  celculeticus  for  the  spherical 
shell  pcteutial  and  the  6-lB  potential  are  in  eacellent  sgreemoat  tilth 
the  erporiaental  data  end  era  oonelderably  better  tbao  either  of  the 
ourvBS  freo  tbs  Oerthelot  equation  or  Piti 


IS  found  to  be  436^^  and  R 


IS  of  de  Rocco  and 


Chlorine  (Cljl 

Free  previous  calcoiatlons,  i 
5,00  A for  Clg,  7he  value  of  r^  was  found  to  be 
poraaeters  Bith  tbe  spherioal  shell  pcteotipl  ti 
Eoover  led  to  the  celculetad  curve  lo  Figure  10.  In  the  ebsenoe  cf 
reliable  second  vlriel  oooffiolent  datn,  second  tirlsl  ooeffioieut 
curves  fran  Pltaer's  oorreletlcn  and  the  Bertbelot  equation  nere  eal- 
oulated  and  are  presented  lo  Figure  10  for  cesnporison.  Tbe  agreement 
at  tbe  Bople  point  is  quite  good  for  all  tbrea  curvea.  Thia  la  rasork- 
oble  uhen  one  considers  the  range  in  (£  valuea  that  Bare  used 

to  celoulate  S /b  for  Clg.  Tbe  ogreHeni  at  loB  teuiporeturee  Is  not  as 
good,  but  is  still  aooapteblo.  If  tbe  Bertbelot  equation  and  Pluer's 
correlation  give  the  correct  vlrUl  ooefficlehto,  then  it  eppears  that 


CBlouIatiofl 


rrcu  UBse  calculations  it  appears  t&at  tbis  eetboi  of 
potactial  anargy  functioas  is  valid,  and  tbat  at  laast  the  sise  parm- 
aters  used  era  satisfeetory.  Use  of  the  aaergy  poraaeters  In  caloula- 
tlcaa  for  Dore  aemplex  Bolocules  should  indicate  ahether  they  are 


Ethane  (CjEj) 

taialag  eethyl  (CB3I  groups,  it  is  necessary  to  have  pari 
group.  lo  viee  of  the  success  of  Claney  aud  Vendersllce 


6-12  potential 


Pitaer  (10}  In  his  eorrelatiou  found  the  acontrio  factor  of  aethauo  to 
he  0.013,  whioh  iudicetes  that  the  petoutiai  energy  function  should  bo 
close  to  a 6-12  potential.  Therefore,  it  seas  resionoble  to  treat 


^haue  as  an  aten  obeying  the 
irs  for  a group. 

Cioaey  and  Vanderslice,  i 


6-12  potential  ii 


The  approach  taken  here 
etba.ne  as  though  it  oere 


s probably  nearly  cc 


the  6-12  potentisl  ore  given 


potentifil 


rj[|3  = ‘>•32  A-  ThB  value  of  f uae  taken  u<  0.T70  A,  nnd  tie 
energy  function  was  calculated  ftoa  Equation  6.  Tkeae  caloulationn  gave 


TCs  naloulated  potential  energy  fiinoticn  for  CgHg  in  preeented  ia 
reduced  fora  In  Table  U and  in  Pigure  11.  A fr-16  potential  ivee  found 
to  fit  the  calculated  curve  very  well,  and  scae  points  are  plotted  In 
Figure  11  for  tbls  potential.  Fcou  Figure  3 it  uas  deteralned  that 


In  order  to  evaluate  it  ia  necessary  to  bs 
paraetara  for  HlrscbCelder,  Curtiss, and  Bird  1 
fton  both  gas  viscosity  and  second  virlal  opeffielent 


leters  frem  gas  viscosity  Be 
' the  6-12  potential  are 


potential 


Exparloantal  saoond  virial  coefficient  dau  for  CjHj  Is  cited  by  Din 
C201.  The  second  virial  coefficient  was  also  oaloulated  frcaPltaer's 
correlation  and  frou  the  Bertholot  equation.  These  values  of  tho  second 
virial  coefficient  ate  given  in  Figure  12. 


Poteotis:  EMrgy  Funcrlcio  of  CjHj 


I 


I I f f' 


'1U8J0JJJM5  BUOMS 


le  expsriEDenial 


Eapirlcally  fUtlae  6-li  pountlal  pi 

One  uould  net  expect  tliie  value  to  be  larger  than 

2.  Tho  calculated  potential  energy  functione  Tor  all  other  aole- 
oulee  have  given  B values  that  were  epprcxiBetely  oorrect.  This  value 

4.  The  value  at  E*/li  ohUlued  empirically  was  about  ISO^,  ee  oca- 
pared  with  230^  fret  gas  viscosity  date,  tfethene  has  an  £*/k  value  of 

pletely  satisfactory. 

The  difficulty  in  fitting  puroiaters  to  these  data  probably  occurs 
because  the  data  cover  a temperature  range  of  only  100*’K,  and  there  are 
no  values  of  the  second  virlal  coefficient  near  the  Boyle  point.  Soall 

side  of  the  peraaoters  under  these  ocodltlcns. 

Since  Pitaor's  oorralatlon  gives  second  virlal  coefficients  that 

calculate  the  ssoond  virlal  oceffleient  over  a wide  range  cf  tempera- 
tures. Using  the  value  of  the  Boyle  point  and  one  other  point  frcai 
Piteer*s  oorrelction  with  the  calculated  second  virlal  cceffloient  for 
the  b-lB  potential  gave  B*/k  s 286°i;  and  B*  : 


6-12  potential 


fclloniag  eslculatioos 
Bsrtbelot  eqnotion  to 


detcraloe  tbe  Boyle  teaperaturs. 


Borthelot  eqnetion 


(E*/l!)j.,2  =1^  a21 
Averegltig  tl 


,t  froD  $ae  Tieooslty 


CB  /k)j_i4  = l.M?  CE'/kJ 

I It  ets  prenloucly  youad  tl 

Sir*-"’ 


r for  the  eetbyl  group. 


( <r*/H|a,  1 r 124"K 

in  ohtalaBd.  Therefore,  the  parasetero  ueod  for  f&3  aro 
£‘/t  z 124‘IC 
r*  = 4.32  k 


Uso  of  tho  epberical  shell  noOel  in  these  cslcuUtions  mthor  then 
Che  6-Ih  pocentlel  would  heve  nede  onlf  ehouc  2 percent  difference  in 
the  flnel  value  of  ^/k.  This  certeinlf  is  well  within  the  aceurecy  of 
tbs  above  celculetioos.  The  caloulated  value  of  ^/h  for  CUu  as  ooly 


Vanderalice's  assiu^tloe  tl 


The  secned  vlrial  ooefficient  calonlated  frcm  these  permieters  and 
Tq  I 3,2a  for  tbs  spherleel  shell  potential  Is  sham  In  Figure  12  and  is 
in  reasonable  agreeaent  witb  those  determined  by  other  methods. 


aiAPIES 


CAlCDUIIOre  FOR  POLYAICmC  UJLECDIES 


1b  Ibis  section  potentlel  energy  lunctioos  and  second  vlrlaX 
coefficients  Here  eeloulaced  for  pclyaccoilo  moleoules  ssiag  os  paraa- 
eters  for  the  peripheral  atcos  the  ateole  paraaoters  calculated  troo 
the  dlatoDle  gee  data.  The  equations  of  Eaniann  and  Laabcrt  ware  again 
used  to  calculate  the  potential  energy  functions.  For  all  interactioas 
•becweea  dissiailar  ataas,  the  folloeing  average  paraaeters  were  used. 


where  and  £q  n: 
and  is  the  peri 
The  auhscrlpta  have 


interactions  between  like  oolooulea 
r interaccioDs  between  oolecules  A and  8. 
aiooaings  with  the  r valaes. 

central  silicon  atra  in  81F.  and 


in  i, 


solved  by  trial  and  error.  Initial  attespta  to  use  the  ssne  parameters 
for  carbon  atoms  as  used  for  fluorion  atone  caused  calculated  values  of 
S A for  the  ..oleonles  that  ware  much  too  largo,  and  it  was  concluded 

and  helium.  Likewise,  the  parsaeters  for  silicon  and  sulfur  atoms 
should  be  boundod  by  those  for  argon  and  neon.  Because  of  this 


Uficsrtainty  in  cboosing  paroisters 


o&Qtral 


Coda  nuuher 

<'*’S  or  Si 

> 

(r-)^  = 3.64  A 

( if* ; 120*% 

4[(r*l*r+ : 3-«  A 

( = 3S.4'V 

!d  the  CSj  parsDstere  given 


For  enaaple,  caleolatlon  SiCCil^l^  - 3 uai 

All  potential  energy  funotioni  uore  caleitlatod  using  ai 
otsputor,  and  the  ainieui  io  oaoh  of  the  potential  energy  ci 
located  by  aaklng  caloalatione  et  intervals  of  0.01  A near  i 
Tsble  16  coDtsins  a suomary  of  the  perBaeters  determined  io  these  cal- 
culations. The  calculations  for  each  nolecule  ore  discussed  separately 


resulted  In  oaloulated  secoad  viriel  ooefflclencs  that  Here  evea  close 
to  the  experimental  values.  Even  Hhaa  ^ I't  for  corhon  mas  tapen  as 
the  same  as  that  for  fluorine  ettoa,  the  calculated  e“/1  »as  sc  saall 
that  the  calculated  second  vlrlal  eoofflolonta  Here  only  about  one-hslf 


surprising  because, 
damonstreted  that  m 


The  value 


Ah  interchtiog  result  of  these  oalculetioos  Has  that  cbaoglng  the 
paraaetsrs  for  carbon  changed  the  valuea  of  E*A,  but  had  little  effect 
on  B or  on  the  shape  of  the  potential  energf  function  when  it  Is 
plotted  00  reduced  ooordlnatea.  Figure  IB  shows  a potential  energy 
function  plotted  from  the  calonlations  of  CF^  - 1,  but  the  other  two 
calculeted  potentials  differed  from  it  by  a negligible  aaount.  For 
this  reason,  only  one  set  of  velires  of  R/B*  end  E/E*  >111  adoguately 
represent  all  tha  calculated  potential  energy  functions.  The  oaleuioted 
values  of  E /k  end  B end  the  reduced  poteotiel  eoergy  fueeticn  for 
theae  oalouletlons  are  given  in  Table  17.  Pointa  plotted  on  Figure  13 
for  a T-21  potentlol  nod  for  r = 1.90  for  the  spherloal  shell  potential 
show  that  aithar  oould  be  used  as  an  ndeguato  aaplrioal  repraaantatien 
of  the  calculated  potential  energy  fusetion  of  CP.. 

In  Figure  14,  the  three  typee  of  interoction  that  aske  np  the 
poteotinl  energy  of  CF^  are  shown  for  CF^  -1.  This  figure  ataous  that 
although  the  interactions  between  osrbon  ond  fluorine  atcas  are  saaller 
than  those  between  fluorine  etms,  they  ere  by  no  weens  negligible. 


These  caloulations  show  why  it  Is  possible  to  aaka  eelcolations 
end  B are  evaluated  eupiricelly,  en 


Flgtire  13-  Fftteatial  Snsrg/  PuncUon  of  CFj 


obtained. 


sccnrate  pbtentisi  energy  fanctloa  can  be 
neglects  the  carboo  atoas  entirely,  the  spherical  shell  ocdel  parameter 
Is  1.90  Ian  exact  fltl,  ubieh  aUs  fits  (enpirioally)  the  potential 
energy  fnnction  vhen  carbon  interactlcas  ere  included.  This  explains 


coefficients  using  the  spharleal  shall  model,  and  of  Cloney  an 
allce'e  Cb}  calculations  of  poteetlal  energy  functions  xlthout 

by  Hamann  and  Lxsbert  (3)  are  shomn,  along  with  calculated  second  vlrl 

sbomn  are  the  oalculotod  curves  using  psrsneters  frco  calculations 
CF^  - 1 and  CPg  - i and  the  value  r z 1-90  for  the  spherical  shell 
ffiodol.  &s  can  be  seen,  the  calculation  for  CF^  - 2 fits  the  data  bett 
tbaa  for  CPs  - 1,  although  both  sra  of  the  right  order  of  magnitude, 
appoars  that  the  calculated  value  of  is  slightly  larger  than  the  b« 
fit  to  the  experimental  data  and  that  the  actual  E*/li  should  be  scae- 
Khera  between  those  fros  CP^  - 1 and  CF^  - 2 calculations.  This  is 
borne  out  by  tbe  fact  that  a good  fit  (by  eye]  to  the  experimental  dal 
'k  z 38i“k,  and  r*  = 1.90  for  tbe 


td  using  R c 4. 
1 potential. 


re  rnnge  of  tasperature,  : 


se  to  the  exporimental  date 
: the  Berthelot  eguation  is  in 
Tbe  second  virial  ooeffioiant 


frma  calculation  CP^  - 3 was  got  as  good  as  far  CP^  1 and  CF^  ~ 2 sod 
Is  not  shown  or  discussed  further. 


Tb«  lalue  of  ^ used  far  the  fluorine  otcns  Her  the  dUteoce  frin 
the  center  of  the  C - C bond  to  the  fluorine  stem,  which  was  calculated 
freo  bond  lengths  ond  angles  for  the  aolecole.  The  value  of  ^ for 

fp  = 1.749  A 
?C  = 0.770  A 

Again,  the  shape  af  the  potentlel  energy  functions  frn  the  differ- 
ent calculations  differed  by  s negligible  Enount.  The  ^ and  E*/k 
values  for  eaoh  culculation  and  the  coemon  reduced  potential  energy 
fnnotian  are  given  In  Table  IB.  The  calculated  potential  energy  func- 
tion is  given  in  reduced  font  in  Figure  16  and  can  probably  be  adeipjntely 
represented  by  using  r^  a 1.6  for  the  spherical  shell  potential. 

ElVerlDental  second  virlal  coefficients  for  CjF^  era  cited  by 
Beed  (21)  and  are  shown  In  Figure  17.  'Also  given  in  this  figure  are 
curves  calouleted  froa  Pltaer's  correlation  and  frui  the  Bertbelot 
equation.  Corves  celonlnted  fron  the  paraeters  of  CjF^  - 1 and  2 
using  B 1.60  for  the  spherical  shell  poteotlal  are  elac  shown.  The 
calonlated  curve  for  C2F^  - 1 agraaa  very  well  with  the  Bertbelot  equa- 
tion end  Pltser's  corroletion  and  appronebes  the  axperlaental  date  at 
lower  teaperatures.  The  calculetod  ourve  for  CgF^  - 2 is  In  belter 
agreeeent  with  the  awperiiaenial  data  at  high  tnnparaturea,  hut  dnviatei 
at  lower  tampernturea.  Nose  of  the  oaloulated  second  virlal  ccefficlant 
curves  predicts  such  a rapid  change  In  second  virlal  coeffiolenta  is  is 
shown  by  the  experLnentel  data.  Although  adjustment  of  the  paraseters 
for  carbon  atesus  could  mshe  the  ooloulatsd  curves  pass  through  the 
enperimental  date,  it  doea  not  appear  likely  that  any  calculated  curve 


fermnns  for  CjFj  for  tbe  7-28  potentUl  enslMttd  fron  gss  ris- 
oooity  Eoasuroaiens*  by  MeCoobrey  wO  Slngfl  (22J  are 

TBU  value  of  E /i  agraes  fairly  wall  with  the  value  852.2'^  fron 
CgF^  - 2,  but  tbii  value  of  E ia  oonsidarably  saallar  then  the  valuea 
5.81  A fren  OjFj  - 1 and  5. 825  A freo  CjFj  - 2,  Aa  might  be  eupeoted 
beeauaa  of  thii  low  value  of  2 , these  paremeters  from  gas  viscosity 
do  not  fit  the  second  vlrlal  ooaffleloot  ee  well  as  oelculotions 
CjFj  - I end  2. 


)t  ehaogo  tli 


The  value  of  ^ used  for  tha  ehlorina  al 
length  1.766  A,  Again,  ohanglng  carbon  pt 
shape  of  the  potential  enorgy  function, 
end  E /i  for  the  throe  oalculations  and  the  cmaon  reduced  potential 
energy  function  are  given  In  Table  Id.  Tha  potontlal  energy  function 
froa  these  oalculations  Is  given  in  Figure  18,  where  polots  plotted  for 
r^  z 1.8  for  the  spherloel  shell  poteneiol  ehcos  tbet  r,,  = 1.6  gives  a 
satiifaotory  fit  ..  the  oelouleted  ourve. 

Ewpeilcental  second  viriel  coeffloients  of  CCl.  have  been  reviewed 
by  Frenoie  and  loGlnshan  (231,  end  tha  ourve  representing  the  enperl- 

Flgure  19,  the  ewperimentel  second  vlrisl  ooefficient  ourve  la  shown, 
along  with  onrves  calonlatad  from  Pitmnr's  correlation  nnd  the  Bertbelot 
agnation.  Also  shown  ere  ourvee  oelouleted  fron  tl 
culatlons  CClg  - 1 end  2 


Figure  19,  Seaoad  Virlal  Coefficient 


Dftntal  s«»nd  virial  coeffieienta  are  available  aver  anoh  a aai 
af  tceiparatura  tbat  aae  ihauld  nat  pUoa  tea  much  iapartence  ai 


Tbe  value  at  ^ tog  tbo  gru^a  was  the  C - C bond  Uagth, 

1.S4I  A.  Tbe  aalculatad  R ead  E /k  for  tbe  ebree  oalculatiaos  ere 
given  la  Table  20.  elong  with  tbe  ctatoa  reduaed  patentiel  energy  funa- 
tlan.  The  aelauleted  pateatial  energy  function  la  sbowa  lu  Figure  20, 
where  it  la  shotm  that  the  apberical  abell  patectlal  far  r„  r 2.10  flta 
the  oalaulated  curve  very  well. 

Expcriuental  deta  an  the  leaand  virial  aaefficlent  of  ClCKj]^  is 
given  by  Raacn  and  Laaibert  (31.  These  data  are  shawa  In  Flgnrn  21, 
along  with  curvea  aaloulated  fren  Pitaar's  carralatlM  aad  the  param- 
otere  fron  aaletlatlon  ClCSj)^  - 1 and  the  aphertnal  abell  potential  for 
r^  = 2.10.  Tbe  Berthelot  eguation  gave  a curve  eaaentially  the  aaae  aa 
piuee'e  eorreletion  end  la  nut  ahonn.  Tbe  agreeaent  aooag  all  curves 
sbCKn  In  Figure  21  is  excellent,  especially  when  considering  that  a 
temperature  ruage  nf  300'^  is  covered. 


CjfCBjlj 

The  value  of  ^ used  for  the  groupa  waa  calculated  frnn  bond 


lengths  end  nogles  end  nss  1.939  A.  The  eelus  of  ^ uniU  for  the  inner 
cerbon  etms  was  0.770  A.  The  oalculsted  e"  and  E*/k  values  for  the 
three  ealculallons  nrd  the  ncsonon  reduced  potential  enerp?  function  nre 
given  In  Table  21.  The  potential  eoerg?  function  is  plotted  in  figure 
22,  along  with  points  for  the  spherical  shell  potential  for  r*  n 1.6 


good  fit  would  prohahly  ha  r.  c 1,65. 


In  the  ahsBooe  of  experiaental  second  virial  ci 
curves  calculated  from  Pltner's  oorralaticn  and  the  Bertholat  eqnotioo 
are  cospared  with  the  curve  for  oalcnlation  - i and  2 parM- 

eters  and  = 1.63  for  the  cpberioal  shell  poteotial  as  shown  in 
Figure  23.  All  of  these  curves  are  in  reasonahle  agreement. 


The  value  of  ^ used  far  the  fluorine  atoos  was  Che  S - P bond 
length,  1.564  A.  Since  ooly  ealculation  SF^  - 3 gnve  roasoueble  results, 
the  potential  energy  function  in  Figure  24  and  Table  22  applies  only  to 
this  CBloulation.  Also  shown  in  Figure  24  are  points  for  the  spherical 


for  the  epherieal  soell  potential  should  fit  Che  calculated  curve  eoiis- 
factorily. 

Enperiaentel  leccnd  vlrlel  ooefflcients  frmu  several  sources  are 
cited  by  Bead  (21)  end  these  points  sre  plotted  in  Figure  25,  The  cor- 
reletloh  of  Pltaer  end  the  Berthelot  egustioh  gsve  preotlcally  Identical 
reaults,  so  mly  the  results  of  the  oelouletlon  fron  Pitaer's  correla- 
tion ore  ahcom  In  this  figure.  Also  given  In  Figure  25  is  the  curve 
for  oalculatlon  SFj  - 3 using  r“  . 1.7S  for  the  sphorlco!  shell 


71 


psteatlQl.  The  oelaulstios  fm  Plizar's 
to  tho  sxpetimaotol  data,  aad  tba  oalotilstod  oarva  far  S'.  - 3 diffara 
trim  tbo  «i;ariBancal  date  by  only  about  10  percent. 

31noe  tha  parraetart  uaod  for  oolfur  in  the  calculation  that  gara 
tha  only  reaaoncble  raaults  u'ero  thoae  for  noon,  sot  thoae  for  orpon, 
thia  ahould  be  dienuaaad  further.  In  all  of  the  molecules  discussed  up 
to  Sj,  the  central  atom  or  atoms  contained  eight  shared  electrons  in 
the  outer  shell.  For  S.,  hohaver,  there  ere  twalre  elastrons  In  tho 


electrons  are  bonnd  more  closely  to  the  fluorine  atoms  ratber  than  to 


The  value  ol 


reduced  potential 

.11  potential  for 


cpiita  large,  changing 
the  reduced  potential 


y fbnotiens  in  Table  23.  The  reduced  potential 
in  Figures  26  snd  27.  The  potential  aaergy 

: 1.0  ss  shoKB  in  Figure  26.  The  potential 


Bafirg?  'luietion  for  Sif . - 2 fits  the  ipberical  sbell  patoatlal  for 
r„  z 1.8  as  sbran  In  Flgata  27. 


Eiperlaeiital  neasuraaiwts  of  the  second  elrlal  oosfficiont  of  Sif^ 

Figure  28,  aloog  sritb  oaloulatod  curves  fren  Plteer's  oorrelatioa  and 
tbe  Bertbelot  equation.  Also  shoNn  are  the  calculated  curves  for 
SiF^  - 1 and  2 using  tbe  values  of  r^  for  tbe  spbarleal  shell  potential 
that  are  given  alove.  Fitter’s  ccrralation  fits  the  data  very  »ell,  aad 


the  Bertbelot  aquetloa  less  vjoll.  The  oalouleted  curves  froa  51F.  - 1 
and  2 were  not  very  close  to  the  etperliiencal  data,  but  one  gives  high 
valnes  and  one  gives  lo»  values  of  the  second  virial  ooefficlent.  These 


tee  eutvee  sppear  to  have  the  sane  sf.ape  as  for  Fitter's  correlation  and 
the  wperinental  data,  uhiob  indicates  that  the  shape  of  the  potential 
energy  function  and  the  velues  of  a"  are  probably  nearly  correct.  There 
la  a very  large  differeace  in  the  E /b  vnluee  fron  the  two  onleulttloos, 
sihloh  is  the  mason  for  tbe  large  diffarenoe  between  the  curves  for 
SIF^  - 1 and  2.  A value  of  E*/k  of  shout  SD0%  would  be  noro  rcasonahle 
than  either  of  the  values  244. 2°K  end  2S2.1^  that  ware  calculated. 

This  Indioetas  that  the  proper  paraneters  for  silicon  in  SIF^  should  be 
alightly  senllar  than  those  for  ergoo,  bat  not  as  anell  as  those  used 


SKCBg),, 

Tbe  value  of  ^ used  in  these  calculations  was  the  Si  - C bond 
length,  1,909  A.  The  cnlculated  peiaatisl  energy  function  for 
Si(CHj)^  - 1 is  given  in  Table  24  and  in  Figure  29.  The  spherical  shell 
potential  for  r*  = 2.90  Is  shown  In  Figure  29  to  fit  the  oalouleted 
potentinl  very  wall. 


fisure  29.  Potential  Ernsr;?  rnuction  ot  Sl(CBj)4-l 


Tbe  sacMd  viri&I  eoeff: 

of  fislcDlatloas  cS  second  vli 
and  tlu  Sertholct  equation  «i 


eat  of  SUCRj)^  haa 
3 and  is  eto«n  in  Figure  30,  the  results 
I oceffielnats  from  Pltser's  corrslatioo 
very  closa  together,  so  oalp  the  results 


of  the  caloulstlca  frm  Pluer’s  correlation  ure  shiwn  In  Figure  30, 

Alto  ihcei  In  this  figure  is  the  curve  celeuleted  free;  the  poreueiers  of 
oelculstics  slfCHgls  - 1 using  r^  = 1.90  for  the  spberloal  shell  oodel, 
this  latter  curve  is  in  onutually  poor  sgreanent  elth  the  experieontal 
date;  the  ealcuUtnd  values  being  about  bO  percent  of  the  siee  of  the 
KtpnriBental  values  of  the  seoond  virlal  coofflotent-  In  vie*  of  the 


calonlatloos 


mPTHt  V 


Siapun  H(D  DISCQSSION 

The  purpose  of  tbe  oalculatious  in  Part  I bos  been  to  oslcuiete 
potenlisi  energp  functions  for  seversi  polyotcmio  molecules  end  tbec  to 

cosffiolsats  end  cooperiog  «itb  experimental  date.  Siooe  It  is  easy  to 
lose  sight  of  the  purpose  of  each  step  in  the  oaleulnttooe  and  how  it 


frta  ioert  gases,  reduced  poteotial  energy  functions  for  the  dlatooio 
gases  Sg*  O2,  ^^2%  ^hre  oelouloted. 

2,  Oeing  these  oalculated  potential  energy  functions  end  6-12 
potential  parameters  frca  experimental  data  on  these  diatomic  gases, 
energy  parameters  for  hydrogec,  fluorine,  chlorloe  atcms  and  methyl 
groups  aero  detsrmined.  The  siae  psrmaeters  chosen  for  the  etoea  io 
these  calculations  were  teeted  by  celculatiag  socand  virial  ooefflcienta 
and  comparing  with  experimental  dale, 

3,  Csing  the  eguatiooe  of  Bmssnn  end  Lambert  and  tha  atcmic  pera- 
eters  determined  from  the  diatmala  gases,  potential  energy  funotions 
were  calculated  for  CF^,  CgFj,  CCl^,  CKHjlj,  CgiCRjlj,  SF^,  SiP,,  end 
SUCH3I,,, 

4,  These  letter  potentiel  energy  functions  waro  tested  by  oeloula- 
tiog  secoad  virial  ccefficleots  and  amapering  with  experlmontal  date. 


Ilieee  oUiMitmiiiia  naTe  mida  by  aaplrioally  titting  sghsilcal  abell 
potmliU  lb  tbt  slleuUted  poteiitlal  ensrgy  fubbtiocs  sad  Ibsa  DsSag 
lbs  tables  of  de  Roooo  and  Hoorer  (6)  to  osloulece  seoond  virial 


Sesules  of  ibass  oslculaiScns  Hill  now  be  dUcusaad  in  aore  detail. 

Calculations  for  tbs  diatomic  gases  sbeaed  that  their  potential 
energy  functions  are  indeed  not  very  Dueh  different  in  shape  from  the 
6-12  petantial.  but  this  differenoo  increases  vith  tbe  sire  of  the  oola- 
cule.  For  exenpLe,  the  potentiel  energy  funotloa  of  hydrogen  could  be 
adegnateiy  represented  by  a 6-14  potential,  while  that  of  chlorine 
reguired  a 6-20  potential.  The  oelculeted  a*  vnlaes  for  the  dlatoaic 


ch  indicates  that  the  sine  of 
different  frm  the  sine  af  tbe 


for  the  fluoTlae  and 
dnta  gave  good  results 


when  used  in  more  ocmpien  molecules,  uhich  indicates  that  the  treatnBBt 
of  dlatcnio  gases  as  dumbell-shaped.  freely-rotating  molecnles  is 


adeguete. 

vara  interaotions  in  the  constltuant  hydrogen  end  carbon  atems  were 

function.  This  result  confirmed  the  findings  of  other  authors  that 

to  a 6-12  poteetial.  Tbe  calculations  for  ethane  end  other  Boleoulee 
containing  methyl  groups  showsd  thst  treating  mathyl  groups  ss  though 
they  were  single  6-12  oenters  of  intercotion  was  an  edegusts  method  of 
approach,  Whi^  egein  eonfirms  the  findings  for  methane. 


Tbi  oabculations 


ocber  polfsttau  aoleeulea  eoatalulas  oarboit 
10  Bjdreseo  «eot)  ebcued  tbit  onij  tta  C - C 


iipe  of  tbs  rsiultlog  poitotiil 
1 did  hoYB  8 larger  offset  upoo 
osleulatlons  Indicates  that  oca* 

) as  in  the  spherical  shell  sodel 
:bat  has  approxlaately  the  correct 
ire  thso  determined  cmplrioally, 
should  he  ohtaiaad. 


In  these  celculatloQS 


obtained  usU;  neon  pamaetera  for  earboo,  although  In  aoas  cases  better 
results  Here  obtained  using  paraoeters  ealdrray  bstween  those  of  neoo  and 
helium.  In  no  oases  aero  helium  psroeters  found  to  be  better  to  use 


for  earboo.  The  parameters  determiaad  frai  the  dlatecilo  gases  mere 
used  for  the  peripheral  atoms,  and  in  rlam  of  tbs  success  of  these  oal- 
eulatlons,  they  must  have  beeu  approxlmntely  eorrecrt. 

There  mas  mora  dlffioulty  In  getting  good  results  for  the  molecules 
oohtaiolhg  slUeen  and  sulfur  as  the  central  etoas  because  the  Intarao- 
lleos  Involviog  the  central  stcns  form  a large  pert  of  the  total  Intomo- 

ouly  ba  obtained  by  using  parsmeiers  that  mere  saiallet  then  those  for 
ergon  (pridubly  ebcut  20  pnrcsnt  smeller),  mhile  for  SF,  It  mas  nec- 
essary to  use  ueoo  parameters  for  sulfur.  This  latter  assumption  mas 
Justified  by  assuming  that  because  of  the  twelve  eleotroos  in  the  outer 
shell,  these  eleotro.ie  ooulti  be  considered  aa  "helanglng”  to  the  fluortne 


8tCQs,  iBsving  nsBaBially  b 
uere  not  siioo«BSfl]l 


tester.  Calonletions  for 
ant  £ /b  psrmeterB 


In  son  OBsse  the  sBlculsted  second  sirinl  CMfriclBsts  did  not 
en&etly  natch  tbs  experiosntnl  data,  but  for  alitost  sll  of  tbose  cnees 
BTUCb  better  egreenent  cculd  have  bean  attained  by  arbitrarily  adjusting 
the  parmeters  eailfsated  for  the  central  atciie.  It  was  felt  that  the 

energy  fusctlono  acre  obtelsed  by  using  a sat  calculation  procedDre.  and 
farther  arbitrary  adjustments  uould  serve  no  useful  purpose  io  tbls 


As  an  interestlog  sldellgbt  of  these  caUuletlons,  it  wei  shorn 
that  Pltaer's  oorrelation  yielded  very  good  estimates  of  the  second 
virisl  coefficients  of  these  scleoules-  The  Bertholct  eguatlon  was 
found  to  give  good  estimstes  of  the  second  virlal  cnefficient  st  Low 
tanpereiutes,  but  was  not  as  good  at  tamperatures  trtmud  the  loyle  point. 

Since  these  calculations  gave  reaooaabXe  rosalts  in  all  but  one 
cssa,  it  eau  he  assumed  that  the  interactive  betneeu  many  simple  poly- 


gualitative 


CKSPJsa  n 


lEASnaElSHT  or  ISl/STIy  AI  tflK  leSIPSSAHJIiES 


The  tbsrail  pc&££nre  coafridien,  OP/dT)^^  hss  baeu  measured 
(29,30,31,32)  tor  several  organic  liquids  at  temperatures  above  The 


la  a pressure  bomb.  7he  liquid  sample  is 
sealed  ia  uith  mercury,  7he  temperature  is 
ura  outside  tbo  mercury  seal  is  adjusted  to 


repeated 


alo  vclune  of  liquid  iaslde  tbe  cell.  Ibis  procedure  Is 
several  temperatures  vitbout  removiag  the  sample,  so  that  a 


series  at  temperatore  snd  preasura  retdtaqs  at  essaatielly  coastaat 


volume  is  Ateised.  Ia  geaeral.  It  mas  fcucd  that  the  pressure  varies 
liaaarly  tdtb  temperature  to  give  ao  observed  value  of  (AP/^TJ,  This 
observed  value  is  then  dorreoted  for  smtli  chenges  la  cell  volume 


arlaiag  free  tbe  tamperaeura  and  presaure  ebaages. 


roaa  temperatura,  so  that  the  temperature  cculd  be  controlled  easily  by 
BOBvunilonel  Essns.  la  tbe  present  «or&,  tbe  thermal  pressure  cemffi- 

B9.S°K,  sad  its  normsl  bolLiag  paint  ia  ec  143.2°!!,  hetb  of  nbich  are 


belOM  tba  freasing  point  of  maroury  aod  well 
Beoaase  of  tba  lor  teEperatura  at  etblcb  tbe  a 
the  present  work,  it  waa  ocoeisary  to  ebaege 


sad  tha  aealicd  of  tesiparotars  oootrol  frect  those  used  by  presiois  expsr- 
Iceaters.  The  use  of  so  eeseacUll;  eossttnt-voluse  cell  v.es  retsiasd, 


fa  la  Figure  31.  Tie  ss 


aicfcol  :ube  ruled  nith  3/16-iaeh  o.d.  heUces  of  24  BI33  aioSel 
sad  Bed  a oapscU;  of  eiiou;  6 oc.  Ike  aiekel  helices  eere  sdded  to 


isctple  tabs  BOS  silver-soldered  to  a aiokel 
brass  beUws-sesl  valve  aas  sersBed.  1 
1 ia  so  eagle  patten,  cod  it  uaa  oecassary  to 
perait  turnios  tbe  valve  rtaotely  frsai  abeve. 


coaaeotios 


saaple  coll.  Tbe  bottoa  of  tbe  sscple  cell  >ss  closed  by  e flexible 
aiekel  ballous.  Tlere  »oro  boOi  Intor.ial  sad  oxieraa!  bollov.s  stops  to 
prevent  deuge  to  ete  bellces  by  over  exteaslao  or  ocr.pression.  A con- 
tact was  so  placad  tbst  the  bellcws  oeda  ocataot  at  zero  pressare  dif- 
ferantial.  Tbe  cootoct  was  conaected  to  so'ladioator, described  by 
Hllla  (331  abicb  was  core  aeasitiva  tboa  tbe  coaveatioaal  electric  laap 
iadioatcr.  A beavytszll  copper  oyliader  eaclosed  tbe  s&nple  cell  to 
fora  a pressure  cell.  Hollua  was  used  to  provide  pressure,  aad  aoterod 
tbe  cell  tbrouab  s 3/!4-ina(i  o-d,  sttl..lees  steel  tube.  Taspersture 
was  aeasured  on  tbe  outside  of  tbe  ssaple  call  by  a tbereocouple  wbleb 
aatered  tbe  pressare  cell  tbrcugb  the  pressure  tubing.  Tbe  wires  to  tbe 
oootaot  also  eatarod  tbraugb  tbls  tubs.  Ccnneotloo  to  the  saaple  eel! 
freet  e glass  vacuuc  systuc  was  else  cade  tbrougb  3/I6-inob  o.d.  staio- 
less  steel  tubing.  T.be  apparatus  was  eaclosed  by  a baavy  brooxe  caa. 


Uvge,  silvered  donar  flaalc. 


Tlia  lid  to  tte  can  aod  all  threaded  boonaetioas  ware  aaaled  by  a lov- 
eelilag  solder  ableb  expaode  upca  frcesriog  to  aote  a tlgbe  oo&nootioa. 

OporatiOJi  of  tbe  apporatua  uas  as  tollna. 

1.  Tbs  ptasstre  sell  nos  evaooated  and  filled  wilt  dry  bslium, 

3.  Liquid  aiirogoa  was  poorad  into  ite  da^ar  flask  so  ttei  it  coo* 
taotad  tte  aetel  oat,  oooling  the  apparatus.  The  apprcocLaate  tenpere- 
ture  of  tbs  can  nea  aaasurad  by  s thsrmooouple  oo  the  side  of  tte  cea. 


traperaturs  and  cooling  bad  slowed  to  8 very  slew  rets,  a saapls  ef 
liquid  was  dUtillad  into  tba  aeapls  cell  froa  a glass  vaceuxi  eystsa  to 


3.  After  ell  liquid  aitrogsn  bed  ei'eporatsd.  aod  the  taaperetare 
begau  to  drift  slowly  upward,  tbe  velve  nes  abut  by  use  of  a rasovable 
rod  ineerted  Into  the  vaive*turtiog  Deshetiaa, 

6.  Frca  50  to  100  pel.  of  bsllLt  uas  adclttad  to  tbe  preeaure  cell 
fros  a oyllnder  of  bigb-presaure  bellut  tbrbugb  a eyetec  of  needle 
valves  wbicb  ooald  be  used  either  to  lowoduoa  heliin  into  or  exbauet 
beiiisB  froo  tbe  pressure  oell.  Ibe  Ittroduetiee  of  beliia  ocspreesed 

7.  Tbe  tespereture  was  allcesd  to  drift  alculy  upward,  end  whsu 
tbe  indieator  signeled  tbee  ibe  bellows  nee  at  the  oouiact  point,  tbe 


B.  Fifty  to  100  pei.  acre  pressure  nos  pul  oa  tbs  celL  end  tbe 


pei.  apparatus  uas  duifraed  to  vltlutantf  1,500  psl,  la  tho  pross:.re 
ooU.)  At  :ni»  tl2e,  more  U^id  nltrosam  ea«  added  to  tbe  denar  flask 
and  the  apparatus  he{;eB  to  oool  again.  Sellom  mas  relaaead  free  the 
Ball  as  tho  apparatus  oooled,  so  that  the  ballCKs  ranined  apprcocIa8tal7 

10.  ffhea  the  apparatus  reached  the  temperature  at  Mhloh  the  valre 


shut,  the  valve  naa  epeaed 
apparatus  or  left  ia  for 


gage  nhicb  mas  factory  calibrated  to  read  ocourately  to  ^,1S  perceot 
of  full  scale  (1,500  pai.l,  i.o.  ±2.25  pal.  gaga  pressure.  The  gage 
pressure  seals  could  be  rood  to  althln  1 pal , , and  preasure  differences 
could  probablf  bemoosured  more  accurately  than  the  ±2.25  pel,  guaras- 


Thermocouple  e.a.f.  values  mere  meaiured  on  e Leads  and  horthrup 

type  K2  potentiometer  upon  mhlcb  e.m.r.  neesuresents  of  ±0.0001  07. 

but  the  rarlation  freen  day  to  day  la  more  than  this  cacunt. 

TemperalucB  oessureoenss  mera  mode  ualcg  eaUiraiod  oopper-ceustao- 
tam  thazmoeonplas.  Xharmcoouple  U-1  nas  inside  the  pressure  cell  on 


matlon  point  of  carbon  dlomide 
oalihrstion  eguatlcc  applies  to 


of  liquid  nitrogen  and  at  the  subli- 
1 ntm,  pressure.  Therefore,  e 
th  thermocouples. 


Tha  c.ethiiiQ  an^  el 
osreehicaUy.  er' 


.breteti  egeiiifib  tbe  vapor  praaeure  of 
ir-presaure  apperarua  of  Crcsidor  CM). 

I Ihe  calibrarioo  were  analyaed  cbrcoab- 
toad  to  oontaie  oolp  traoo  Ijapurftlaa, 


sample  purity  aboald  bavo  caaaoo  oai[il 

1 vapor  pressure  points  for  tbsraoooupl 
t z - W.&67  + 0. 231890  - O.PlSOio^ 

I t z Taaparsturo,  9c 


Axte-jion 


of  tbe  ceil  bec« 


nou^reprodccib^e  ccnbbQb  indicator. 

U cftraaialy  it«gblgibia.  7ha  apting 
eaougb  ca  tbt^  graaauro  .rrora  becaa 
also  bo  aaslaobod. 

Corrections  ioT  Ibe  ob&ngo  In  a 
aura  aoo  prassuro  ehaagas  can  be  r.s 
Tbo  obange  la  Intanal  volnne  of  tba  cell  la 
&7  = VO^aT  - Vp^&F 
vnaro  z Change  in  volujse  of  call 

CCn  = Coafflclenc  of  tbernal  er<;onUea  of  aocal  in  cell 
jSi  ■■  CmprottIbSllty  of  -eta!  da  cell 
If  s 0.  than  T cbservod  would  be  oorreot  end  equal  to 
cap/Sl)y.  Since  tils  is  not  true,  extra  pressure  AP  - APj^,, 
bnra  bo  be  applied  to  the  liquid  to  occpensese  for  AV.  Thet  i: 


«=7.-vp 

r Cbservod  value  of 
^ r Ccapressibillty  oi 


Wien  Equations  3S  aud  SC  arc  ocnbiaad, 


Calodlations  showed  t..ac  the  uouor  of  liquid  laaids  the  brass 


-8oo“c-=t«s;6%  ua) 

‘ ! + 66t  H-  0t2  + I'-.S  + #-.« 

where  the  eeitecaats  for  nielcel  ere 
ei=  12.5*  10-4  »C-l 
^ = 8,75  X 10-« 


S:  6. 25  X 10-15 

Tl;!9  «qnatioa  xst  utnl  to  oaloiiUto  «!1  ecbiool  oosffiolenU  of  tiiotatl 


Tbe  aost  otltloal  footor  is 
K&s  the  fssipsratur*  cootrol.  Si: 
ihBsa  oxoeriaoots  do;«aded  upoo  t 

aesta  for  ssch  of  tcaoe  cospousds  sill  be  dlocusaed  separstol^. 


Toohaiosl  gsEiie  CXl^  t:ae  distillod  ir 
Stine  ooloaa  equisalent  to  sixty  thsoretic 
o-laptasa  : sothyloyelGbexaufi.  X oonlwr  o 
ttisse  expatioests.  Tbe  vapor  prassuvs  of 
with  a cslibratad  therfflocoupla  sod  was  fou 
vapor  prosaurs  data  for  CCl^.  Vfiis  isdlia 
ououeh  for  ass. 

tegporatsro  rlsa  was  0.01S°i;/zin.  ia  eis  a 


I to  opree  with  publishad 


sbflpld  bsvs  boos  aasosiiall 
rapidly  with  toipersurre,  s. 


arezBsi  asd  0.013°C/alo. 
‘0  slow  anough  so  that  tiie  apparatus 
It  aquillbriiSB  throsgboui  t] 
la  aids  of  tba  valve  differed  fi 
.ojt  0.4«C,  bat 
Its.  Sisoe  (SP/SI)^  doei 
1 tccperetura  grtdioots  ti 


Tbe  pisBsura  teapcraziira  values  fav  the  two  ceasureueuta 
ata  givaa  la  Table  Cl  of  Appendix  C.  Tbe  praaaure  «as  found  lo  vary 
lioaarlf  iritn  poaparzluze,  and  strazobi  liuea  were  fitted  to  tbe  data  by 
tbs  xatbvd  of  least  iquavae.  Tte  values  of  flP/AT  observed  Here 
obtained  free  tba  slopes  of  the  Hues,  std  uKe  toxpevsturaa  at  ehl^  tbe 


Tbere  Is  aeae  variation  io  the  cboloe  of  the 


teaperature  at  ehlch  a 
Teaperatutes  that  have 


1.  The  vapor  prcsscre  of  tha  subsUaoe. 

2.  The  tmparaturo  at  uhioh  the  ebsolnte  preoscre  la  aaro. 

2,  The  teapereture  at  nhloh  the  cace  pressure  is  zero. 

The  last  xethod  xas  used  To  the  present  xarX  beoeusa  of  tha  ease  In  cce- 
putatlOB,  The  difference  it  the  oethods  is  puite  soail  and  cao  be 
aegleeted  for  vapor  pressurea  up  to  the  rorxal  boiling  point  of  tha 

The  oorreotioa  to  be  edded  to  the  .eeasured  values  of 

(9P/0T>^  nes  found  tn  bo  S.42  eta./°C.  Tha  results  of  these  experl- 


Thls  uorx,  Bennlnga  & Soott  CsO} 

(9P/ST)v,  OP/eTly  Porcent 

Aa,/‘C  Ata./“C  Doviatlnn 


These  resulte  are  In  eueellent  ecreeaent  aith  those  of  Bonninpe  end 


k ffuoh  flora  Ifl^rtast  tvaporit' 


aade  eofiCaot. 

jr«  srodlont  oxlsteti  Slone  Isngth 


sturo  ef  tla  volrs  nss  £rmi  2 to  3*%  Malar  tba-n  tSa  teaparalurs 
flsasurea  inrlds  the  pressure  eell,  It  ras  sot  poesiMs  to  dotamlos  the 
eredioots  aloag  the  sflflple  oell  laclde  the  prossnre  cell  heosose  oaly 
one  therfflooojple  oould  be  loaielled  is  the  pressure  cell.  Since  the 


pressure  cell  ess  etde  c£  copper,  nliioh  hos  high  theroal  oonductl/-ty, 
the  tecpsretura  elvoa  the  lenatli  of  the  scspls  oell  should  have  been 
fairly  soifero.  If  flcst  cl  thie  prsdlect  toot  place  la  tha  valve,  than 


ij  ebwt  5 iwreant,  tr>  wan  If  thote  m 

tba  teaparctura  v;aa  eaaauraf.  Sicca  xbe  tt 
Isctely  at  the  center  of  tiL«  actpla  call,  a 
racdlcB  ahculd  Cara  beaa  oCtalcaC, 
prassuca  caafficient  aboiUd  Ca  buoC  battar 
coaid  icdicBte.  TMa  gradlact  aas  aocatacl 
tCat  oCa&sas  lo  gradiaata 

Conpacy,  and  sac  9S  parcact  pure  or  bac 

Hally  pura.  ISa  aajor  Icparity  «B8  CO 
in  an  axpariaactal  aeaauraaact  a.iaiad  o 
after  pcrfficatica. 


sal  praascra  oi 
i a Sdir  taoparacura  gcadicat 


;Can  CCa  gradleot  of  3% 
during  tCo  ceaauraaienta,  ao 
aupariitantal  arror. 
puroCaaad  iroa  tCe  Hatnaaon 
aa  auppliad.  Tea  caterial 

ind  analyalc  of  aatarial  uaad 
a traoa  of  CO2  racalcad 


licaarly  tH:C  (.Co  a.a.f.  of  tba.moca.'pla  ri-1  khlsC  war  at  tea  oactar  of 
tea  aeaplB  call.  TCa  Gtaarured  praaaure  nas  ploccad  agaiaat  tCa  a.a.f. 
of  tCaisoooupla  b'-l  aad  tie  alopa  red  Iniarcopt  at  aaro  page  praraure 
i»ra  datarcicad  grapCically.  Tlaaa  alcpaa  and  intercepta  cotMoad  with 

oSaarrai  raluas  of  is/i?  and  teaparituro,  TCa  axptrlreotal  valuer  of 
a.a.f.  and  proaacre  ara  tabulatad  it  Appendix  C. 

TCa  cherral  prasauro  ooaffioiaot  of  a anbatasoa  it  ralotcd  to  ixa 
eonpraasibility  a:d  tie  cccfficfanc  of  tCaraal  expaaaion  by 

=s£/&  as 


Tlie  value  of  c4  coo  be  deterained 


Ima  lltiali  densitr  data.  Id  tbe  next 
DolDr  voliUDtt  of  liquids  ran  be  corre* 


Proa  liquid  deneity  dura  fo 
ooustants  in  this  equation 


B z 2.0159  X 10-® 

so.*  X “K 

Tbeee  values  of  Xcan  be  used  uith  tbe  observed  values  of  ^P/aT  to 
estlaate  ^ , stbich  is  Deeded  to  taebe  tbe  correctioua  te  .&P/AT  for 
thersal  expanaiou  of  tbe  cell.  After  corrected  values  el  (8P/9Tly 
have  bean  obtolnod,  tbey  osu  be  used  witb  tbe  values  of  ct  to  cslculate 

culatioDS  are  given  in  Table  25  and  in  Pigures  52  end  53. 

lu  Figure  32,  tbe  tbenaal  pressure  eoefficieet  of  CF^  is  shown,  and 
tbe  curve  tbroug]:  the  points  was  dresvo  by  oyo  as  being  tbe  best  represen- 
tation of  tba  dote.  There  are  three  points  that  deviate  noticoebly  frea 


Figure  32.  Jkerosl  Pressure  Coofflelent 


PIfiure  33,  CnprasslbUlt?  of  CP. 


The  reaalti  of  tDeae  neasureaent*  wiil  ba  fiiacuMed  further  iu 
later  cbaptera. 


Direct  experlmeatal  ueasureaeats  of  the  theraal  prcasure  oeeffi- 
worfc-  The  ratea  of  toopareture  riae,  taEporewre  groUlents,  and  pre- 

Phillips  Research  Grade  fflsthaao  «es  used  for  these  experiseote  and 
uas  99.60  pereeat  pure  as  supplied.  Traces  of  air  and  oarhoa  dioxide 

The  results  of  these  neasnteaaats  of  (5P/9T)e  for  CR.  are  gisea 
io  Table  26  and  in  Figure  34.  The  rax  exporlueatal  data  are  given  io 
Appendix  C. 


Thd  liquid  ctaprufisibility  can  be  caloulatsd  from  Telocity  of 
soucd,  best  oepeoity,  end  deiieitj  data  for  a liquid.  CooMninq  tie 
liquid  Mnpreeaiblllty  aod  tie  coefflolent  of  tberoal  expanaioc  yield* 
tba  cberoal  presaare  coeffioieat.  Sovllnaoo  (3d),  la  his  tables,  dees 
values  of  the  theraal  pressure  coefficient  celeulated  in  this  naoaer, 

oonsiderable  difference  beteeoo  the  direct  Deesurasents  of  this  work 


this  difference  mounts  to  about  IS  percent,  because  of  tbs  cueber  of 
different  quantities  to  be  determined  for  the  Indirect  method,  the 
direct  neesurenents  in  the  present  work  ere  preferred. 


was  evcntirally  found  to  be  a defect  in  the  coavoiutioaa  of  thn  nickel 
belims.  The  leak  rate  woa  about  7 x 10*^  moles  of  helina  per  hoar 
with  1 atm.  pressure  difference  across  the  belltfss.  This  leak  rate 
probably  would  be  saaller  with  other  gases.  Since  the  pressure  dif- 
ference ecrcss  the  bellows  Is  always  small  during  operetlco  of  tba  appa- 
ratus. there  should  be  liitlo  or  no  effect  by  this  leek.  The  bellows 
were  In  e cesapreased  condition  most  of  thn  time  during  e oeaeurment  so 
that  if  there  were  any  net  flow,  hellun  would  be  forced  into  the  cell. 

If  haliicn  were  forced  into  the  bell,  it  would  either  I'orm  e bubble  or 
go  into  solution.  If  a bubble  were  formed,  the  value  of  tSP/SIlv 
would  be  that  of  s ges  oed  would  be  much  lower  then  the  values  observed. 
Helium  is  not  very  soluble  in  liquids  st  the  tmperatures  used  In  these 


t&s  good  resDlto 


oooorlbutad  aotbiDg  i 


CH4PTER  VII 


ISE  SMOOTHED  PCTEKTIAL  CEU  UOESL  OP  THE  LIQUID  STAR 

Theie  hav«  bean  eeverel  Bbdels  of  tbe  li^iS  itate  that  have  been 
uaed  to  corcoUts  and  calculate  pbyaioal  and  chermodynanic  proportlea  af 
liquids.  Host  of  thesa  aiodela  are  quits  ooaplex  and  can  only  ba  heodlad 
adequately  by  numerical  oeeaa.  The  end  results  of  these  labers  hove  not 
been  espeoiaUy  gretifylng,  end  because  of  the  difficulty  in  using  these 
ccaplcm  models,  this  author  does  not  consider  thee  to  be  appropriate  for 


Therefore,  the  model  that  mill  be  applied  to  the  liquid  state  iu 

tlODS.  Calculations  in  this  chapter  mill  shorn  the  nodal  to  be  quall- 

the  other  liquid-state  models,  the  quantitative  egreauent  eitb  real 
liquids  is  not  good. 

The  call  model  of  the  liquid  state  uses  the  assumption  that  each 

oell  that  is  on  the  order  of  the  molecular  sine.  The  force  field  ectiog 
upon  each  mcleoule  is  rapidly  fluctuating  beceuse  of  noleeular  motion 
and  can  ba  replaced  by  a field  of  apbarical  eymetry.  The  simple  cell 


aiBlBBuIes  axcbange  placBB  with  a Irequaaoy  that  laoraaaaa  with  tha 
vhluae  par  sileciila.  Tiiua,  the  alBpla  cell  woilel  oacnn  ha  uiah  far 
calculatloa  of  triBSport  propertiea.  If  the  average  leegtb  of  tiaa  a 
aolecale  spaods  la  a cell  ta  icffUlentlp  long,  this  chaining  of  placet 
will  have  little  effect  upon  the  thetaodynanio  properties.  The  cell 
Dodel  givei  wore  order  to  the  liquid  stale  then  la  eotually  precent,  but 
thia  error  is  not  aa  Important  at  1(«  tamperatoraa  and  high  deneity  aa 

llquida  in  the  region  of  the  oritical  point  in  the  preient  work.  There 
are  aevaral  raraiocs  of  the  ceil  model  ivhloh  are  diaouaaed  by  Kirsch- 
feldar,  Curtiss,  and  Bird  (1),  and  by  Frlgogina  t2).  Further  disonaaion 
in  the  present  work  will  be  limited  to  the  smoothed  potential  cell  model, 
as  presented  by  Prlgogloo. 

The  smoothed  potential  oell  model  la  based  an  the  assumption  that 
the  average  potential  energy  of  isteraotion  of  e molecule  in  a cell  can 
be  replaced  by  the  value  of  tbe  pntectiel  energy  when  nil  moleeulea  are 
at  the  oentera  of  their  cella.  The  potential  energy  of  the  molecule  at 
all  poaitlona  within  the  oell  la  aaaumod  to  be  equal  to  thia  average 
value,  ncd  tbe  molecule  la  not  permitted  to  move  outaide  the  cell.  In 
order  to  preaeot  these  asaumptlona  more  clearly,  It  is  necessary  to 
define  them  mathematically.  The  symbola  used  are 

wlr)  S Potential  energy  of  Uternotlon  of  a molecule  with  its  neighbors 
Cewpressed  es  energy  per  molecule) 

w(D)  s Potential  energy  wben  the  moleools  is  nt  the  center  of  its  cnll 


nstants  whleb  tslce  non-aearest  aalgDbor  iDtereotions 
QOn-fltareit  neighbor  intoractioofi  are  nagleoted. 


conatanu  la  the  faoa-canterati 


lattioa  etruoinra  for 
ad  to  deteroilne  theae 


6-12  potaotlal.  tbaaa  conataou  ere  <L,2) 


potential, 


Kolar  Volima  of  LiqoUs 
Ac  vaniabing  extern 


further  oalculationa  involving  thaae 


the  etpiation  of  atnte  for  a liqalO 


(Equatloa  44) 


folloving. 


TliereforB.  frcoi  this  e^uatloa  of  state  the  7o1ubb  per  ooleoiile  Bt  abso- 

tlon  ia  equal  to  S . Ooly  wbBo  noo-neBreet  naiglibaT  iatereotioas  are 
neglecteO  does  the  yolume  at  absolute  aero  equal  r . 

Bead  (211  found  tbat  the  taore  oouplex  1000800-^0008,  DeroastLire 
equation  of  state  for  the  6-12  poieatiel  predicted  tbet  a plot  of  1/V^ 

liquid.!  Equation  45  oannot  bo  solved  explicitly  for  i/v^,  and  ft  oau- 
not  be  seen  by  Idspootlon  that  the  equatioo  of  state  fron  tba  smoothed 
potential  oell  model  predlota  this.  One  can,  boNOvar,  calculate  numeri- 
cBlly  the  values  n/aE  that  oorresponds  to  given  v/v  values  and  plot 
these  values  to  determine  whether  this  Dodal  predicts  that  l/V^  T 
should  plot  linearly.  This  was  dons  for  the  6-12  and  T-25  potentials 
and  the  results  ore  ebon  in  Figure  35.  In  this  figure  it  ean  bo  seen 
that  the  plot  is  exactly  straight  for  the  6-12  potential  and  that  it 
curves  slightly  for  the  7^6  potential.  The  curvature  for  the  7-28 
potential  Is  so  slight  that  it  would  be  undateetable  over  the  range  of 
temperatures  usually  oewerod  by  most  liquids  up  to  their  normal  boiling 
poiata.  Thus,  according  to  this  modal  the  plots  of  lA^  against  T 
should  be  very  nearly  linear  for  molecules  that  have  very  narrow 


Plg-ra  35.  S^uatio-i  of  State  it  Zero  Extant:  P.*3suce 
far  ttsa  ^oattai]  Pateatlal  Cell  t'cCel 


ctfiDtial  Bcerg)  fuDcti 

Xa  Figures  2b  aad  37.  1/ 
veluM  e{  tlia  taalar  volune. 
Area  tba  refarances  ^Ivaa  ba3 


18  as  io  tha  7-2S  potantiol  as  sell  as  for 
^ is  platted  against  T using  oxpsrlaental 


fz 

CCCHj)^ 


turs  St  high  teaperaturas  where  the  aasimptions  of  the  cell  model  are 
DOC  as  good  as  at  lower  temperatures.  Values  of  the  eanatasts  A and  B 
in  Eguetion  40  were  evaluated  for  the  cempounds  abeam  in  Figures  3b  and 
37.  These  values,  as  well  ss  the  extrapolated  molar  volumes  et  ahsoluie 


liquids  at  absolute  zeri 
methods.  For  CH.,  Kurt: 


0)  are  in  verv  good  agreament  with  the  valu< 
ir  hydrocarbons  using  other  eKtrspoletiOD 
eve  the  value  23.0  cc./mole  which  oemphres 


Id  I-C^i2,  Kurt*  gi 


le  value  40.33  froa  the 
S4.1  00. /mole  while  tl 


molecule  iu 


[•aa/siQLj  'jOt 


poteatlal 


ir  vibratioas  at  absoluta  zara  ivill  eeuee  zona  addliiaa  to 
tbe  TOlima  at  abaoluta  zero  over  wbat  ia  ahom  la  Equation  it,  Evaa  if 
tbeaa  Tibratiooa  ara  not  lakan  Into  aooount.  tba  value  of  v*/v  at  0°K 
la  only  1-092  for  the  6-12  potential,  wblcb  aaJcaa  R /d  = 1.03  at  0*^. 
For  tbe  T-28  potential  at  abaoluta  zero.  Equation  <6  qivea  v*/v  z 1.031 
and  R /d  5 1. 01.  Tbarafora,  the  aolecular  aeparatim  at  abaoluta  zero 
abnild  giro  eatimatoa  of  R that  are  lo»  by  a nazlBum  of  3 percent.  For 
oloae  packing  of  spbarioal  taoleoulaa  at  abaoluta  aero,  tba  Dtolacular 
aaparatioa  ia  (1,3} 


oeaparlaon.  It  oio  be  aeon  Ibet,  in  general,  the  valuea  of  S*  hm 
Equation  47  are  raallar  than  both  tbo  6-12  potential  ualuaa  fraa  oxparl- 
oantal  data  and  the  oaUulatad  valuea  froa  Part  I.  The  calculated 
veluea  fron  Part  1 aaro  generully  loeer  than  the  6-12  voluaa  acd  ara  lo 
better  agraament  with  tbe  a"  valaos  froa  Equation  47.  Thaae  reaulta 


Batter  egreenent  Kith  the  R°  valuea  frm  tha  aotar  volune  data 
pcidtably  oonld  have  bean  obtained  in  aoue  caaea  by  ampliioally  fitting 
to  ozperiaental  aocond  vlrlol  coefficient  data  8*  and  e"  volaes  for  the 


calculated  aUapcs  oi  the  potential  energy  funotiona  from  Pert  I.  Ifhea 
thia  eaa  done  icr  CP.,  tbo  veiue  B*  z 4.70  A eaa  obtained,  whloh  ia  in 
better  agreement  vlth  the  value  4.S4  A [rm  nolar  volunte  date. 


Thermal  Preaaure  Coefficleot  .of  Liguide 

The  cbeiaal  preaaure  ooeffinieot  for  the  smoothed  potential  cell 
aodel  cee  be  determined  by  partial  differentiation  of  Equatioo  44  to  be 


IPT  V = 


..  1 


It  is  inteceuting  to  note  frcm  this  equatioo  that  the  only  dependenea  ol 
{dP/9T)y  upon  the  e>tponentB  o and  n in  an  o-n  potential  oooura  in  the 
paraaecer  a,  mhicb  is  not  stroogly  seneltlve  to  changei  in  m and  a. 
Boarrangcoent  of  this  equation  ahaua  that  a plot  of  ^V(9P/87)yj''^ 
against  should  be  a stroight  line  with  elope  Nha  and 

intercept  1/Kk.  Note  thot  V /e^  repraaenta  the  molar  volume  of  the 
liquid  uhen  the  average  separetion  of  the  molooulea  is  equal  to  Bg.  Th4 
value  of  fv”l^^a  shaild  be  leas  than  the  cube  root  of  the  molar  voliai 
at  absolute  sere,  but  of  the  same  order  of  magnitude,  (Por  the  6-12 
potential  s s 1,1224,  and  for  the  7-2B  poteotieL  a b 1.0683.1 


to  Pigures  28,  29,  Bad  40,  tbaimal  preaaure  coefficient  data  for 
A,  Og.  Ng,  CCl.,  CB4.  and  CP4  are  plotted  as  described  above.  The  data 


for  A,  1^.  Og.  0014,  I 


Rcmlinson  (38).  The  other  volues  of  OP/8T)y  for  CH^  and  the  values 
for  CP4  are  the  values  determined  enperlmentnlly  lo  the  present  norh. 
Thare  is  seme  curvature  in  the  lines  dram  through  tOe  points  in  these 
figures,  but  the  slopes  do  not  change  rapidly..  Since  there  is  sme 


ft]  the  figures, 


jie 


deteraioed  fron  the«e  plots.  It  is  laterestlog,  boMover,  to  see  1£  tbs 
slopes  of  the  carves  are  of  the  order  of  oageltude  that  would  be 
expected  if  the  voluaes  of  the  supercooled  liquids  et  absolute  xero 
were  used  in  place  of  V /s^.  The  voloes  correspoeding  to  V*,  V®,  aod 
V /a  sre  in  the  crder  V*/s^^V®«v".  Straight  liaes  with  slopes 

are  shown  for  ceaparison  in  these  figures.  (These  slopes 
ore  slightly  larger  in  absolute  value  than  they  would  be  if  t’/s^  had 
been  used.)  The  intercepts  of  the  actual  curves  were  larger  than  1/Kk, 
and  the  dashed  lines  of  coastant  slope  ware  arbitrarily  positioned  on 


the  ordlnato  for  easy  cooipariaoa  of  slopes. 

The  values  of  the  thermal  pressure  ooaffioient  of  argon  in  ICMlln- 
son's  tables  is  a oonbinacion  of  direot  maasurewents  above  «0°K  and 
values  below  90%  oalculated  frcm  velocity  of  aound  data.  Those  values 


oulated  valuei  hi 


values  oslculsted  frcss  the  velocity  of  sound  at  lower  taperatures. 

This  Dshes  the  value  of  the  thermal  pressure  ooaffioient  change  rapidly 
with  teaparature  at  the  lower  temparatuiaa  in  the  tnbies  of  Benlinson. 
The  slope  of  the  curve  in  Figure  30  for  the  direot  nensurements  in  in 
such  bettor  agreement  with  the  elope  then  for  the  curve 

teken  frot  Rowlinson'e  tables. 


Hitt  this  explsnation  of  the  curve  for  ergon,  the  slopes  of  the 
curves  fron  axperimencal  end  calculated  thermal  pressure  coeffiolenu 
given  in  Figures  30,  39,  sod  40  are  in  good  egceaaent  with  tha  slopes 
taken  from  the  liquid  noler  volume  at  abaolute  ware  tenparetnre.  The 
curve  from  the  axparimental  measurmnenu  of  the  presant  work  for  CO.  la 


as  exceptioa  to  tMo  ogreeaent.  In  vIdu  of  tho  good  agrooDeat  between 
the  ectuBl  slopes  and  the  slopes  fron  for  all  the  other  conpouncts 

it  sppeers  that  adOitlonal  Beetureneats  for  CR^  Bight  be  edrlsahle. 

B for  Liquids 

Asother  quantity  that  can  be  oaloulated  froe  the  smoothed  poteotlal 
cell  Bodel  la  x,  whiob  is  defined  as 


By  differentiating  Equation  43  for  the  snootbed  potential  cell  model,  x 


the  7-^a  potentials  and  these  values  were  plotted  to  give  the  curves  in 
Figure  41.  These  curves  show  that  over  tbe'range  of  v*/v  normally 
occnrring  In  the  liquid  i0.d«v*/v<0.91  state,  the  value  of  x should  be 
larger  for  anbstanoes  whiob  follow  the  7-28  potential  than  for  those 
thot  follow  the  8-12  potential.  It  is  probably  true  that  x becmea 
larger  at  a given  value  of  v /v  as  the  reduced  potential  energy  function 


becacos  narreerer  (as  □ and  n inereasel. 

The  value  of  x for  real  substacees  can  bo  determined  from  Equation 


oonfigurational  energy  data  is  detertaloed  froa  vapor  pressure  data.  To 
oaloulate  CdEg/dVlj  would  snount  to  differentiating  the  experiooutel 


It  sfieas  reasaaa&le  to  aooueo  that  ooo  be  used  to  opproxlaats  V 
it  deterRloing  reduced  liquid  voluaes.  lo  Figure  d2,  x is  plotted 
sgeioet  VA^  lor  tbe  liquids  meationod  ebooe. 

Tbe  aneothed  poteatlal  cell  model  predicts  that  x sbould  iaorease 

aeoeseerili  so.  Since  values  of  x include  moesureoents  of  E^,  V,  and 
( dP/STly  Mben  the  thermal  pressure  coofficieats  ere  measured  direct  1/ 

ealcuUted  iadirectZy,  x can  be  determiaed  ao  better  thee  to  within  1 
peroeat.  Therefore,  too  much  Importance  should  aot  be  given  to  the 
shapes  of  the  curves  la  Figure  42. 


potential  oell  model  that  at  the  seme  value 
narroner  reduced  poteatlal  eaergf  functions 
of  X.  The  potential  eaergr  functions  of  CF 


prediction  fraa  tba  smoothed 
of  v/v  , molecules  with 
should  have  higher  values 


on  reducod  coordlaates  then  those  for  the  other  matorials  shown  in 
Figure  42,  Hildebrend  end  Soott  (411  give  values  of  x for  other  nonpolar 
polyacealo  molecules  which  would  bn  expected  to  have  fairly  narrow 
reduced  poteatinl  energy  funotiona.  Some  of  these  veloea  ere 


n-CqBjg 

Stars 

SnClg 


These  values,  being  relatively  high. 


0 CE^^,  Sx^riAAatsl  Csta,  sl:U  worEc 

^ Cil^,  Fm  RottliEESoa  (S3) 

Cxporlcextsi  Dstc,  this  vork 


Fl{TJre  t2.  x as  a Fasation  of  Lf^olC  folms 


I 


coaolUBiCRS  given 


regarding  the  shape  of  the  potential 


The  thaoreticel  ilgnifioancB  of  x is  sooeohat  aohalous.  The  egue- 
tion  of  stnte  of  van  der  Baals  prodleta  that  x should  be  unity,  so  thet 

0 von  der  Keels'  fluid.  Since  the  total  range  of  x values  Mas  froe 
about  0.8  to  about  1.1,  van  der  Keels'  aquation  of  state  should  ocae 
fairly  close  to  fitting  oonflgnrational  energies  for  the  liquid  state. 
The  oeleulatioas  in  this  ohepter  have  shown  that  the  saoothed 
potential  cell  csodel  gives  good  qualitative  predictions  about  the 


about  the  potentiei  energy  functions  of  the  aoleculas,  bowever.  Par 
eXEsple,  it  was  found  thet  the  molar  volumes  of  liquids  cen  be  correla- 


ted by  the  aeae  type  of  eqvetion,  regardless  of  the  potential  energy 
functioa,  and  the  same  Has  found  to  be  true  for  the  thensal  pressure 
coefficient.  The  value  of  x,  bouever,  did  give  sane  indioetion  of  the 
width  of  the  redcced  potential  energy  fnnotian,  and  thn  predictioa  of 
the  smoothed  potentisl  cell  model  that  x should  be  larger  for  moleoulea 


their  reduced  potential 


funotinns 


caA?TEa 


CWrEHIEillONAl  ENBEK-VOLCE  EEUTIONSBIP  IN  U90WS 


The  distribution  of  noopoler  oflleoules  about  « given  nolecule  in  a 
pure  ll(piid  Kill,  m the  auarage,  be  spherioally  syaietrloal.  If  kb 
select  I differential  volime  element  dr  at  a distance  R froa  the  given 
molecule,  then  the  probability  of  finding  a molecule  in  tbia  volume 
elemsot  is  tlygiSldr,  Kbere  n„  is  the  average  number  of  moleoules  pec 
unit  volume,  and  gOl)  is  the  radial  distributimi  function  Kblob  describes 
the  density  variations  in  Che  liquid.  Actually,  the  radial  distribution 
function  is  also  dependent  upon  the  molar  volume  end  the  temperature,  so 
that  it  is  proper  to  »rlto  me  rodial  diairlbutlon  fanctinn  as  gtB,V,T). 

In  terms  of  ibe  radial  distribution  function,  the  oonfigurntionui 


If  one  assumes  that  the  mutual  pair  potential  enargy  function  is  of  the 
Lennnrd*voeos  m-n  type,  then  actually  UCB3  is  a fuootioa  of  E*  and  R/B* 


os  a genorel  form  of  this  type  of  potential.  Also,  It  has  frequently 
been  nnsiaed  tbst  g(B,7.T)  depends  only  upon  the  potential  energy 


energy  puraiaeteri  for 


that  or.e  cokld  tbsn  Ju$t  be  ejbII  write 

jcs.v.T)  = ay 


Kaking  tbese  EulMtltijtl«i<  le  Equation  Si,  ml  ehaagiug  the  variehle  oC 
iDtegretien  front  E to  y gives 


i«2eV3 


2nJ  Oj(y)Bi(y,^.i*lAr 


patentlsl  energy  fonotioo  of  the  nclecuies,  the  aolar  volisoe  along  the 
sBiuration  curva  of  a liquid  can  be  expressed  la  terms  of  T*.  The 
rlgbc-hond  side  of  Equation  57  is  then  actually  a fuaotion  only  of  T° 
along  the  saturation  curve.  The  left-hand  side  of  Equation  57  is  the 
negative  of  the  reduced  van  der  Keels'  a s -a°.  From  Equation  57,  one 
can  detereino  that  for  molecules  which  have'  the  same  form  of  poteotial 
energy  fuaotion  (same  values  of  a and  a},  plots  of  a*  against  T*  shculd 
give  a single  curva.  0ns  would  not  uecessarily  expect  that  a*  would 
be  the  some  for  molecules  which  have  poteotial  energy  functions  of 
different  shapes,  however. 

In  the  derivation  of  Hildebrand  and  Hoed  Ml, 43)  which  is  used  in 
developing  Bildabrand's  theory  of  solutitms,  the  assumption  is  made  that 
a is  the  SBse  for  all  suhstancoa.  Since  this  is  an  important  point  in 
the  theory,  it  Is  desirable  to  see  if  a is  really  the  seme  for  all 


nfigurativisl 


calBiUstad  frca 
tloaal  eaarglas 
daacribsd  balw. 


B of  Rwlinson  fd8}.  CtIculotloDS  of  conflgura- 


CP^.  An  equation  for  tbe  vapor  preaeuro  of  CP^  is  oited  b;  Bead 
<21).  and  hast  of  vaporlzatioo  valuea  Mate  oaUulatad  froa  tbia  aqua- 
tion uaing  the  Clapayroa  equation.  Tha  coafiguratioaal  anerpy  ana  than 
oalculatad  using  Pitzar's  correlation  (10)  to  detarmioa  the  enthalpy 
m the  real  and  ideal  gasas  and  to  oorrect  tha  vopor 
0 liquid  aolar  voluaa  vas  calculated  froa  the  equation  pre- 
he  prerious  chapter. 


CgH^.  The  eathalpy  of  the  ide 
range  end  tha  heat  of  vaporization 
given  by  Din  (20).  Also  given  are 


real  and  idaal-gas  eathelpy  for  the 
given.  These  data  tfare  lufflolent 


aal  gas  over  a portion  of  tha  liquid 

the  liquid  voluuas  over  tha  liquid 
sed  to  oeloulate  tha  dlffarauoe  between 

to  oeloulate  S*. 


CgFfa.  Bead  (21)  oites  tao  values  of  the  enthalpy  of  vaporisation 
boiling  point.  Pilsar's  oorralatioo  ass  used  to  oslculote  the  differeooe 


ClCll^ln.  The  best  of  vaporisation  at  two  teuiparatursa  and  tba 
liquid  oolsr  voliuse  are  given  in  Befereooe  2B,  Pitsor's  oorreiatloo  aaa 
used  tu  caleulato  tha  differeuoe  betseen  reel  aad  idasl-gss  enthalpy. 


v&Iuas  of  liquid  ooufiguriiionol  escrg/  ajid  oolar  volumes 

In  Pisuro  d3,  o (dlBSnsionless)  is  plotted  sgalsst  for  tbe  oos'' 
pounds  mentlonsd  sbovo,  noiug  &-12  potential  parameters  for  E and  R . 
The  psrccetars  used  for  this  figure  are  given  la  Table  30.  Tor  the 
BOlaoulas  A,  N2,  O2,  CU4,  end  C23j,>hloh  eooie  fairly  close  to  fitting 
the  6-lE  potential,  the  value  of  a*  is  reas^iehly  close  to  the  same 
(oboot  a 10  percent  rengel  at  a given  I . for  the  moleoules  CP4,  C^Fg, 
CClj,  and  CICE3I4.  ehloh  mere  showi  in  Part  I to  have  potential  energy 
fuootions  that  Bare  very  cmoh  different  frcm  a 6-lE  potential,  the 
value  of  a*  is  very  much  different. 

Since  all  of  the  EDoleeulea  above,  even  the  diatiaic  gases,  do  not 
folloR  the  d-12  poteotial.  the  Bide  range  in  the  values  of  a*  might  he 
thought  to  occur  heccuse  of  this.  In  Fignre  44,  a*  is  plotted  against 
t”,  but  this  time  the  potential  energy  functioa  pereoeters  used  were 
the  ^ sod  E values,  celouleteJ  in  Part  X.  whloh  gave  the  beet  agree- 
meet  with  axparifflantol  second  virlal  coefficient  data.  The  values  need 
are  given  in  Tahle  30.  Use  of  these  parsoiMers  oeHted  the  vslues  of 
s*  for  the  various  compounds  to  heve  Jost  as  wide  a range  as  when  the 
6-12  paroeiers  were  used.  Although  the  calculated  potential  energy 


molecules  which  have  different  shapes  of  potential  energy  functions. 
Toe  Hildebrand  theory  of  solutioca  abould,  as  suggested  by  Bead  (44), 
be  modified  to  tahe  the  difference  io  the  veluee  of  a*  into  account. 
Since  a*  is  evidently  dependent  upon  the  shnpe  of  the  poteotial  energy 


a 


psrmetor  tbst 


potential 


eoergr  ftuictlon  into  eccouat  Caucb  as  Pitaer's  LS  eight  poaaUily  be 
usea  to  cake  ao  leprovemeat  in  the  BiUebreod  tbeorp. 


CHAPTER 


SUmART  OP  PAST  H 


In  Cfaapter  VI,  oeasureDentn  of  tRn  therasl  pressure  eoefneieor  of 
CCl^,  CF^,  und  CH^  Mere  described  end  tbe  results  Kere  presented.  Tbe 
results  for  CCl^  nproe  with  otber  published  euperloentel  meosurraents. 
bat  tbe  results  Cor  CB^  ere  in  poor  egroeaent  wicb  values  celculated 
froB  velocity  of  sound  dots.  Tbo  euperliseiitsl  resolts  of  tbe  present 


sork  sre  preferred  f( 
couparlscoi  tiltb  tbe  results  in  the  present  work  foi 
Tbe  mootbed  potential  call  eodel  ues  applied 
and  it  das  predicted  that  plots  of  1/V^  (V  n jsolar 
d be  approxleately  linear  for  all 


liquids.  This  wes 


te,  and  fra  plots  of  experlnental  molsr 
tbe  supercooled  llqaid  at  obsolute  sere 
severe!  Odupounds.  These  extrapolated 
used  to  estimate  values  of  R for  tbe  s 


sidered.  Values  nf  ■ oeloulsted  fron  7®  mere  senerslly  smaller  then 
R values  obtained  enpirioally  for  tbe  b-12  potential  and  sssller  than 
tbe  R*  valaes  for  the  aora  resliitle  potentials  calculsted  in  Port  1. 

bettor  ooreement  mltb  ibe  R 's  at  absolute  aero  teeperoture.  Tbe  R*’s 
oslculated  frea  V®  values  sbould  be  smaller  then  R*'s  fron  the  oaloula- 
tlons  nf  Part  I by  op  to  3 peroent,  but  the  sotual  differenoe  is  larger 


Part  I might  be  alightly  larger  t 
With  the  smoothed  potential 


therael  pressure  coafflolent  should  also  be  tnaeositlTe  to  changes  in 
potential  energy  {unctiona,  end  that  plots  of  [vOP/ST)^J~^  against 

end  some  curvature  mss  found.  The  slopes  of  these  curves  Here  found  to 
be  In  good  agreosient  altb  trhet  aould  be  expected  fros  the  smoothed 
potential  cell  model  booed  upon  the  molar  volume  of  the  oompounJs  at 
absolute  aero  temperacure. 

The  ffiiootbed  potaotial  call  model  else  predicted  that  the  value  of 
should  be  lorger  for  molecules  nhlcb  havs  narrmver 


lo  Chaptar  VIII,  it  mas  ahcEsn  that  the  value  of  the  reduced  van  dor 
Uaals*  a 9 a is  not  the  seme  for  ell  moleoules,  bot  dapcnda  upon  the 
shape  of  the  potential  energy  funetien.  This  finding  suggests  that  s 


roexamisatiou  i 


of  the  Kildebr^d  theory  of  aolutiona. 
same  for  all  aubatances. 


Although  the  daleuletioos  lo  Part  II  are  intarasting,  thay  give 
little  guaotilative  iaformation  shout  the  potential  energy  functions  cl 
the  molecules  studied.  The  results  of  Part  I regarding  the  shnpes  of 
the  potential  energy  funations  are  aupported  by  the  results  of  Part  II. 


gualitative  predlotionc 
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SPHERICAL  SHELL  ECIZNIIiU. 


SPHERICAL  SHELL  POTEmAl 


lEWARD-JONSS  6-16  POIEmlAL 


r/r*  0/e'“  r/r'  0/f‘ 


la 


i/r*  v/e° 


164 


D/f° 


fSSSSCSE  CCESnCSSKT 


IKEaKM. 


192 


e CCeFTICIEKI  OP 


193 


;S8sg 


1 


W Cf4 


:§^§B8S3S 
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PsriReter  defined  by  Elation  1 
Geme  function 

Value  of  ~U(r)  at  tbe  niniaia  In  Q(r).  Usually  used  as  a 
p Distance  froo  the  center  of  a oeleeule  to  the  center  of  a 


in  Pltaer'a  cerrelation  (1 
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